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the Navy Operational Regional Atmospheric Prediction System (NORAPS). 
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slightly less than a quarter wavelength is most favorable. The initial meridional 
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shear vorticitv, which depend on the structure of the mean How, are important factors 
in determining the most favorable cross-track tilt. In the jet streak case with strong 
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An east coast cvclogenesis event during 19-21 January 1986 is used to illustrate 
the effects of superposition. Despite strong upper-level forcing, rapid cyclogenesis at 
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vorticity center and the surface cyclone was not achieved. Modifications to the initial 
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I. INTRODUCTION 



Maritime extratropical cyclogenesis events are among the most powerful and 
destructive weather systems outside the tropics. The western oceans are prime areas 
for cyclogenesis, especially during the cold season (Miller, 1946; U.S. Navy, 1955; many 
others). Because most of these systems intensify steadily, they are relatively well- 
predicted. However, about 20 percent of the cyclogenesis events include a short 
(12-36 h) period of extremely rapid intensification. This intensification, which is often 
poorly forecast, results in a rapid rise in wind velocity and sea heights. Improved 
understanding of these "rapid deepeners" will be beneficial to coastal forecasters, the 
commercial shipping and fishing industries, and to the U.S. Navy. 

For tactical and strategic reasons, it is often necessary or desirable for naval task 
groups to operate in regions where winter storms are common. Accurate and timely 
forecasts of the onset and intensity of these events are required so that operational 
commanders can' make decisions regarding ship and aircraft deployments and select 
appropriate missions. Although hazards to mariners from high winds and seas are of 
primary concern, many decisions are influenced by direct and indirect effects of 
cyclones. Important tactical variables that are greatly influenced by the passage of a 
strong storm include radar or electro- optical sensor performance and the acoustic 
characteristics of the upper ocean. 

Unfortunately, the storms that are most dangerous and have the greatest impact 
on naval operations are the least well understood and the least reliably forecast. 
Numerical weather prediction models have a poor track record in predicting the 
magnitude and timing of intensification in these rapidly deepening systems (Sanders 
and Gyakum, 1980; Bosart, 1981). To better forecast these lows, it is necessary to 
achieve a better understanding of the interactions between the synoptic-scale, 
mesoscale and boundary layer processes that are required for rapid intensification. 
Sanders and Gyakum (1980) arbitrarily chose a deepening rate of 1 mb h' 1 for 24 h 
(adjusted for latitude by the factor sin (p/sin 60°) as the division between ordinary 
cvclogenesis and rapid deepening. At 40° N, this is approximately equivalent to IS mb 
in 24 h. This deepening rate has been generally accepted as the lower end of the rapid 
deepener scale. The most rapidly intensifying storm during the three-year period 
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studied by Sanders and Gyakum had a deepening rate more than three times the 
minimum value. 

A general hypothesis of rapid cyclogenesis has been defined for the Experiment 
on Rapidly Intensifying Cyclones over the Atlantic (ERICA), which is part of a five- 
year program of basic research sponsored by the Office of Naval Research (Hadlock, 
19S7). In condensed form, the hypothesis is that rapid cyclogenesis in the western 
Atlantic depends on the proper superposition of strong positive vorticity advection 
(PVA) aloft over a pre-existing low-level cyclone in a favorable environment. There 
are several contributing factors that are related to this hypothesis. Some of these 
factors may also be present in cases which are not rapid deepeners. . 

Jet-level processes may enhance cyclogenesis in several ways. Maximum 
divergence aloft, which encourages upward vertical motion and convection, is expected 
when the low-level center is under the left-exit or right-entrance quadrants of a jet 
streak. In the exit region, the lower branch of the ageostrophic circulation assists low- 
level convergence into the surface center and provides a source of moisture for deep 
convection, which encourages coupling of the upper and lower systems and increases 
the baroclinic transfer of eddy available potential energy (EAPE) to eddy kinetic energy 
(EKE). 

Jet streaks also advect positive vorticity into the upper levels above the surface 
system. Even with straight flow at the jet level, an intense jet with strong cyclonic 
shear on the poleward side can provide significant PVA. The magnitude of the PVA 
will be increased even more if the jet streak is propagating through the base of a short 
wave trough or from trough to ridge. 

A pre-existing maritime cyclone provides a starting point for rapid deepening. 
Presence of a pre-existing perturbation reduces the time required for the growth of a 
small-amplitude perturbation to a significant amplitude cyclone. The combination of 
PVA aloft, a surface low, and a favorable environment must be maintained for some 
time if the new system is to become sufficiently organized to achieve rapid deepening. 
The surface low circulation also organizes low-level moisture over a relatively deep 
layer. 

For rapid cyclogenesis, the low-level environment must be highly baroclinic and 
have low static stability. Such an environment is created when a deep layer of cold air 
flows over warm surfaces, especially in regions with large thermal gradients. These 
conditions are most frequently satisfied following winter season cold outbreaks over 
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warm oceans. Strong surface temperature gradients can be found at the land-sea 
interface or over the "north wall" of the Gulf Stream. Cold-air damming along the 
Atlantic seaboard can greatly intensify the low-level baroclinic zone along the coast. 
Synoptic patterns that cause cold outbreaks over the western Atlantic are also 
favorable for subsequent cold-air damming events. The optimal environment for 
explosive deepening in the western Atlantic is hypothesized to be one in which cold-air 
damming intensifies the coastal baroclinic zone, while cold air flows steadily 
southeastward across a strong sea surface temperature gradient at the edge of the Gulf 
Stream. Continued flow of cold air over the warm water will maintain low-level 
baroclinity and lead to large air-sea fluxes of sensible and latent heat. Thus, these 
processes may work together to broaden and strengthen the rapid deepening domain. 

In summary, many factors can play important roles in rapidly intensifying 
cyclones. This study will concentrate on the importance of the "proper superposition" 
between the source of PVA aloft and the existing surface low. For this research, 
rapidly intensifying cyclones ("rapid deepeners") will be defined to have intensification 
rates equal to or greater than the Sanders and Gvakum definition. 

In this study, finite amplitude -perturbations at lower and upper levels will be 

placed in a variety of superpositions in two atmospheric models-a two-level analytical 

* 

model and a channel version of the Navy Operational Regional Atmospheric Prediction 
System (N'ORAPS). The emphasis* is on the development of /m/re-amplitude 
disturbances that are placed in a baroclinically unstable mean flow. The perturbations 
are inserted at 400 mb and 800 mb in both models and the subsequent development at 
both levels is described. As only one perturbation wavelength is present, this approach 
differs from classical experiments in which infinitesimal perturbations with a spectrum 
of wavelengths are allowed to grow in a baroclinically unstable mean flow. In this 
study, the perturbation wavelength represents the most unstable mode that eventually 
would dominate in the classical experiments for the mean conditions in each case. 
Different combinations of perturbation amplitudes and positions are tested to 
demonstrate the influence of superposition on cvclogenesis rates. One immediate goal 
of the study is to provide further evidence in support of the superposition aspects of 
the ERICA hypothesis stated above. A better understanding of the superposition 
aspects will ultimately improve forecasting of these dangerous maritime cyclones. 
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A review of literature pertaining to rapid cyclogenesis is presented in Chapter II. 
The idealized model results are in Chapters III (two-level model) and IV (NORAPS 
channel model). In Chapter V, the conclusions from the model studies are tested in a 
case study from the Genesis of Atlantic Lows Experiment (GALE). Chapter VI is a 
summary of the conclusions from the research. 



II. LITERATURE REVIEW 



A number of papers pertaining to rapidly intensifying cyclones have been 
published since the landmark study of Sanders and Gyakum (1980), which first 
separated rapid deepening from ordinary' cyclogenesis. Some of the relevant works are 
reviewed here, along with some older references that are also applicable to the study of 
rapid cvclogenesis. The references will be addressed as they apply to the key 
requirements of the hypothesis: 

• upper-level PVA; 

• pre-existing low-level cyclone; and 

• favorable environment. 

A. UPPER-LEVEL PROCESSES 

The superposition hypothesis requires that a source of positive vorticity aloft be 
properly positioned with respect to an existing surface cyclone that has a horizontal 
scale that closely matches the scale of maximum baroclinic ihstabilitv of the jet-level 
feature. Both jet streaks and mobile short waves may serve as sources of upper-level 
PVA. 

The Sanders and Gyakum (1980) study included 267 cases of rapid deepening 
that occurred in the Northern Hemisphere during the three cold seasons from autumn 
1976 to spring 1979. The average position of the surface low was approximately 400 n 
mi east-northeast of a mobile 500 mb trough. This mobile upper trough was found in 
the southwest quadrant of the surface cyclone in more than 75 percent of the cases. 
More generally, the nearest 500 mb trough was in the western semi-circle in over 90 
percent of the events. In a study of 48 cases of rapid deepening in the western North 
Atlantic Ocean, Sanders (1986) concluded that rapid cyclogenesis is fundamentally a 
baroclinic process, and that a prominent absolute vorticity (£ a ) maximum at 500 mb 
contributes significantly to rapid deepening. 

This scenario of a propagating upper trough overtaking a deepening surface 
cyclone is consistent with the classical Petterssen development equation. Petterssen 
(1956) defined cyclone intensification as the change in geostrophic surface vorticity, 
and related it to the change in upper-level vorticity and the change in thickness 
between the upper-level pressure surface and 1000 mb. Thus, it is expected that 
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surface cyclogenesis is associated with the region of positive vorticity advection (PVA) 
in advance of the mobile 500 mb trough. 

Cyclonic shear vorticity on the poleward side of the jet enhances PVA above 
perturbations on the cold side of the surface baroclinic zone. This may explain Sanders 
and Gyakum's (1980) observation that the maximum in explosive cyclogenesis tends to 
occur just north of the strongest tropospheric baroclinic zone. Based on data from 
nine winters, Blackmon et al. (1977) also found that the band of maximum baroclinic 
eddy activity was below the cyclonic shear zone on the poleward side of the strongest 
jet streams. 

In addition to upper-level short waves, alongstream variability in the intensity of 
the jet is important in cyclogenesis. The polar front jet is characterized by propagating 
cores with higher speeds. These isotach maxima were defined by Palmen and Newton 
(1969) as jet streaks. The jet streaks propagate eastward much more slowly than the 
winds at jet level. Therefore, parcels transiting the streak accelerate in the entrance 
region and then decelerate in the exit region. The thermally direct circulation in the 
entrance region, with divergence aloft in the right quadrant and convergence in the left, 
converts available potential energy to kinetic energy. As the parcel exits the streak, 
kinetic' energy is converted to potential energy by the thermally indirect circulation. In 
the exit region, there is divergence aloft in the left quadrant and convergence in the 
right. 

Uccellini and Johnson (1979) demonstrated the importance of jet streaks in the 
development of severe convection. The upper tropospheric jet streak advects cool, dry 
air into the upper part of the cyclone environment, enhances upper-level divergence, 
and transports sensible heat dow'nstream from the convective region. Additionally, the 
low'er branch of the ageostrophic transverse circulation below' the exit region of the 
upper jet streak can induce a low-level jet (LLJ) that imports warm, moist air into the 
convective region. This low-level jet, termed the "conveyor belt" by Harrold (1973), 
can provide a connection between a source of w T arm, moist air (such as the w r arm 
Atlantic Ocean) and the convective region. The combined result of the upper jet streak 
and the LLJ is to create a region of convective instability with low-level convergence 
and upper-level divergence. 

Shapiro (1983) suggested that proper alignment of the upper and lower jets can 
facilitate deep convection and couple the low'-level system to the jet-level processes. As 
an upper-level jet overtakes a surface front, the pattern showm in Fig. 2.1a is formed. 
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The LLJ advects warm, moist air into the area, while the upper jet provides cooler, 
drier air aloft, and makes the lapse rate convectively unstable. The ascent branches of 
the upper and lower secondary circulations are vertically aligned (Fig. 2.1b). which 
produces a deep convective plume just ahead of the surface cold front. 




Figure 2.1. Vertically coupled upper and lower jet/front systems and associated 
secondary circulations, (a) Jet streak exit situated above the surface front and 
low-level jet. (b) Cross-section along the line BET of (a). (Shapiro, 1983) 



Calland (19S3) studied a rapid cyclogenesis event that occurred in the western 
North Pacific during 12-16 January 1979. Because cyclogenesis occurred under almost 
straight 500 mb flow, this was a case in which upper-level support for development (in 
terms of PVA) was expected to be minimal. However, vorticity budget analyses 
showed that horizontal vorticity advection was important, even without a significant 
short wave trough aloft. A large eddy vorticity flux at upper levels was associated with 
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advection of cyclonic shear vorticity by an intense jet streak. Rapid intensification 
occurred when the surface low moved under the upper divergence associated with the 
left-exit quadrant. 

An Atlantic polar low that also developed under relatively straight flow at 500 
mb was studied by Cook (1983). The incipient low formed south of Iceland on 26 
January 1979 and deepened rapidly as it moved southeastward. As in Calland's case, 
upper-level eddy vorticity transport into the budget volume was present despite the 
lack of an appreciable short wave aloft. The rapid deepening coincided with the 
movement of the divergent left-exit quadrant of a propagating jet streak over the 
surface low. 

The Presidents' Day storm of. 18-20 February 1979 has been extensively studied 
by Bosart (1981), Bosart and Lin. (1984) and Uccellini et al. (1984, 1985). The surface 
low formed in an inverted trough along the South Carolina coast and moved 
northward along a coastal front. The cyclone was influenced by two upper-level jet 
streaks. Genesis of the low in the coastal trough was accompanied by an intensifying 
subtropical jet (STJ) streak moving from Kentucky to the Virginia/North Carolina 
coast (Uccellini et al., 1984). A strongly ageostrophic low-level jet over the Carolinas 
advected moisture northward and upward, and contributed to a region of moderate to 
heavy snow in the western Carolinas. Rapid deepening in the Presidents' Day storm 
occurred later, when a polar front jet (PFJ) streak and associated short wave trough 
moved eastward over the coastal low. 

Uccellini (1986) argues that upper-level baroclimc processes also contributed to 
the rapid intensification of the September 1978 QE II storm (so named due to the 
damage suffered by the oceanliner Queen Elizabeth II). Gvakum ( 1983a. b) and Anthes 
et al. (1983) attribute the QE II deepening to diabatic processes, but Uccellini showed 
evidence that an intensifying short wave trough/jet streak system upstream of the 
surface low could have contributed to the rapid development. 

Numerical studies have also emphasized the role of upper-level jets in rapid 
cvclogenesis. Mullen (1979) composited NMC analyses for the genesis periods of 22 
eastern North Pacific polar lows. Hodur (1984) used this composite jet. with maximum 
winds approximately twice the values found in a mean winter jet, to study 
intensification of Pacific polar lows. He found that the linear growth rates associated 
with the mean winter jet were approximately 0.54 d* 1 , while growth rates associated 
with the genesis jet were 1.2 d* 1 . Hayes (1985) suggests that intense jets crossing the 
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coast from continent to sea are conducive to development of surface lows with 
horizontal scales typical of maritime cyclones. Sinclair and Elsberry (1986) describe a 
case of cyclogenesis in a numerical simulation that occurred on the anticyclonic side of 
the jet. Nevertheless, jet-level processes appear to have been important, since the 
development took place under the divergence associated with the right-entrance 
quadrant of a propagating jet streak. 

The importance of jet processes was also highlighted by Ebersole (1984) and Peak 
et al. (1985). Forecasts made from all available data were compared to forecasts made 
without satellite-retrieved information. For the Calland (1983) case, the "no satellite" 
forecasts were surprisingly superior. The difference was traced to discrepancies in the 
position and strength of the upper-level jet stream. This sensitivity to the strength of 
the jet stream reinforces the idea that jet streaks play a crucial role in explosive 
cyclogenesis. 

B. THE PRE-EXISTING CYCLONE 

The second key factor in the hypothesis is the presence of a pre-existing maritime 
cyclone at low levels. In a composite study of rapidly intensifying Atlantic cyclones. 
Rogers and Bosart (1986) showed that rapid deepeners typically evolve from’ shallow 
lower-tropospheric systems. The presence of this low between the Appalachians and 
the Gulf Stream allows rapid cyclogenesis to begin immediately when a baroclinically 
unstable jet streak of proper scale is superposed. 

In the absence of a pre-existing surface low. cvclogenesis will still occur when a 
baroclinically unstable jet streak crosses the coast to a favorable low-level environment. 
This unstable basic state encourages exponential growth of small perturbations, with 
the wavelength of maximum growth being shorter over the ocean than over continents. 
The decrease in the most unstable wavelength is due to reduced friction, absence of 
orography, reduced static stability, and fluxes of sensible and latent heat (Tokioka, 
1973; Gambo, 1976; Staley and Gall, 1977; Sandgathe, 1981). 

This exponential growth from perturbation to incipient low requires a finite 
amount of time. Since the "rapid deepening zone" has a limited east-west extent, the 
favorable superposition of upper support and conducive environment has a limited 
lifetime for each case. Without a pre-existing system, the available time within the 
zone is only sufficient to develop a new finite amplitude low from the small initial 
disturbance. The new low has then already propagated out of the rapid deepening 
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domain before rapid cyclogenesis can occur. Therefore, a pre-existing maritime cyclone 
that is in the proper position is a prerequisite for explosive deepening. 

Two winter synoptic patterns were identified by Miller (1946). Type A cyclones 
form as a wave on a cold front that divides warm, moist air over the ocean from cold, 
dry air associated with a large anticyclone centered over the Midwest. These cyclones 
typically form over the ocean and move northeastward along the existing low-level 
baroclinic zone. 

Miller's Type B cyclogenesis takes place southeast of an older low situated in the 
Midwest or Great Lakes region. The new center forms near the coast, often along the 
warm fro.nt of the old system. A wedge of cold continental air propagates southward 
along the eastern slope of the Appalachians from New England and retards the warm 
front over the coastal plain. This cold-air damming against the Appalachians 
(Richwien, 1980) does not slow the northward movement of the front over the ocean, 
and the front becomes oriented NE-SW along the land-sea boundary. The secondary 
cyclone can then form as a wave on the convoluted warm front to the south, or along 
the coastal front. Kocin and Uccellini (l985a,b) studied a 24-year sample of 18 east 
coast snowstorms. Fifteen of these cases eventually qualified as rapid deepeners. In 
nearly every case, cold-air damming and coastal frontogenesis occurred. 

Orographic processes can also favor cyclogenesis east of the Appalachians. Cold 
air damming associated with the Appalachians is often a precursor of Type B 
cyclogenesis. The mountains can also influence cold frontal passages from west to 
east. The surface front is impeded by the relatively high orography in the Carolinas, 
but proceeds more rapidly across Georgia. This causes a wave-like inflection in the 
low level front (Type A cyclogenesis). In the middle and upper troposphere, the 
frontal zone continues its eastward traverse above the mountains, which results in 
static destabilization that promotes convective activity and enhances cyclogenesis. 

Westerly flow over a mountain results in an induced ridge over upwind slopes 
and a corresponding trough on the downwind side (Holton, 1979; Hayes, 1985). The 
downstream trough can enhance cyclogenesis. A similar effect may be caused by cold- 
air damming. A thermal ridge between the Appalachians and the Atlantic with a 
thermal trough over the warm water would effectively move the area of enhanced 
cyclogenesis eastward. 

Recent studies of western Atlantic explosive cyclogenesis events all discuss rapid 
deepening from a pre-existing cyclone. The Presidents' Day storm of 1979 started as a 



wave on a coastal front over the Carolinas (Bosart, 1981). Of the 15 cases of rapid 
deepening documented by Kocin and Uccellini ( 1985a, b), 12 are classed as Miller 
Type B (secondary low forming on the warm front/coastal front) and the remaining 
three (Miller Type A) moved to a position just off the east coast before rapid 
intensification commenced. Finally, the QE II storm started as a shallow baroclinic 
system over the New Jersey coast (Gyakum, 1983a, b). 

C. THE FAVORABLE ENVIRONMENT 

The third requisite for rapid cyclogenesis is that the superposition between the 
pre-existing surface low and the upper-level jet streak takes place in an environment 
conducive to cyclogenesis. The overwhelming tendency for rapid deepening to occur 
over oceans (Sanders and Gyakum, 1980; Roebber, 1984) provides ample 
circumstantial evidence that the most favorable conditions are maritime. Several 
physical processes have been suggested to explain this correlation. 

During the cold season, the ocean acts as a huge reservoir of heat and moisture 
for the atmosphere. • Cold continental air flowing over the relatively warm water gains 
heat and moisture through upward fluxes of sensible and latent heat. Gall and 
Johnson (1971) calculated the generation of available potential energy by sensible 
heating during North Atlantic cyclogenesis to be at least 34 percent of the boundary 
layer frictional dissipation. Bunker (1976) and Bunker and Worthington (1976) showed 
that the upward air-sea fluxes during winter' are a maximum over the warm water south 
of the Gulf Stream. These sensible heat fluxes are due to the large air-sea temperature 
difference. Additionally, cold season mean winds are higher than warm season values 
(U.S. Navy, 1955), which enhances the heat and moisture fluxes from the ocean. 

The direct effect of the upward energy flux is to destabilize the lowest layers of 
the air. Sensible heat warms the bottom of the air mass and decreases static stability. 
Concurrently, the flux of latent heat increases convective instability. Nuss and Anthes 
(1987) found that fluxes of sensible and latent heat in the warm sector updraft 
contribute significantly to cyclone intensification. 

A number of analytical and numerical studies of general cvclogenesis have 
emphasized the important effect of decreased low-level stability. When the Richardson 
number (Rj) of the mean flow is small (low static stability and high vertical wind 
shear), the growth rate increases and the wavelength of maximum instability decreases 
(Tokioka. 1973; Gambo, 1976; Gall, 1976a; Staley and Gall. 1977; Sandgathe. 1981). 



Low Rj and increased convection induced by air-sea fluxes are also important in the 
genesis of polar lows (Reed, 1979; Mullen, 1979; Sardie and Warner, 1983; Hodur, 
1984; Orlanski, 1984). 

In case studies of rapid cyclogenesis, warming and moistening of the boundary 
layer appear to be important. Both Calland (1983), for a North Pacific case, and Cook 
(1983), for a North Atlantic polar low, noted a dramatic decrease in low-level static 
stability during the early stages of cyclogenesis. In the Presidents' Day case, Bosart 
(1981) attributed the location of the strongest surface winds (in the northern semi- 
circle) to cold air flowing over warm water, which sustained superadiabatic lapse rates 
near the surface. Gyakum (1983a) cites increasingly deep layers of potential instability 
as being important in the evolution of the QE II storm. 

One important consequence of destabilization by air-sea fluxes is a marked 
increase in cumulus convection. The same processes that make the environment 
favorable for rapid deepening (sensible heating of the lowest layer, an unlimited 
moisture source, latent heat flux) also result in convective instability. Thus, it seems 
that 'rapid cvclogenesis and cumulus convection are closely related. Cumulus 
convection has been showm to enhance baroclinic energy conversions from zonal 
available potential energy (ZAPE) to eddy available potential energy (EAPE) and from 
EAPE to eddy kinetic energy (Gall, 1976b; Hayashi and Golder, 1981; Robertson and 
Smith, 1983). 

A number of studies have connected rapid cyclogenesis with the effects of 
convection. Pyke (1965) and Simpson (1969) proposed that increased cumulus 
convection can contribute to low-level convergence, w'hich increases the vorticitv. 
Bosart (1981) stated that a CISK-type process could not be ruled out as a contributing 
factor in the 1979 Presidents' Day storm, wiiile Bosart and Lin (1984) discovered that 
the vertical advection of vorticitv was a key factor in the same storm. Gyakum 
( 1983a, b) suggested that convection was important in the evolution of the QE II storm, 
and that a CISK mechanism was responsible for a notable increase in potential 
vorticitv during the rapid deepening. Anthes et al. (1983) numerically simulated the 
QE II storm, and found that cumulus convection and air-sea fluxes of heat and 
moisture were important in the development. Rogers and Bosart (1986) also concluded 
that cumulus convection enhances intensification. 

Convection may play an important role in forming the ''bridge'' between the pre- 
existing low’ and the upper-level jet. Mullen (1979) and Hodur (1984) found that 
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convection was required to explain the vertical extent of polar lows. However, not all 
cases of rapid cyclogenesis include significant cumulus convection. Mullen (1982) 
documented rapid deepening in a relatively dry air mass over south-central Canada. 
Although the low-stability, high-moisture environment of most rapid deepeners 
encourages cumulus convection, rapid deepening can occur without these conditions. 

Another mechanism that apparently makes the maritime environment conducive 
to rapid cyclogenesis is enhanced baroclinity in the boundary layer. The Richardson 
number, which is used as a stability parameter in several studies mentioned above, is 
also a measure of baroclinity. The vertical wind shear term in the denominator is 
proportional to the horizontal temperature gradient. Consequently, baroclinity in the 
lower troposphere contributes to lower values of the Richardson number. 

Sanders and Gyakum (1980) presented evidence that rapid cyclogenesis occurs 
over the largest gradients in sea-surface temperature (SST). They also suggested that 
larger deepening rates in the Atlantic are related to SST gradients that are double those 
of the Pacific. Sanders (1986) found that the most rapid deepening in the strongest 
Atlantic cyclones occurs as the surface lows cross the extremely strong SST gradient 
along the northern edge of the Gulf Stream. Pandolfo (1985) suggested that the 
generation of cyclonic potential vorticitv at low levels might be significantly enhanced 
by SST gradients across Gulf Stream meanders and eddies. Sandgathe (1981) found 
that locating an upper jet over the zone of largest SST gradient was conducive to 
cyclogenesis. 

Roebber (1984) identified a local maximum in occurrence of rapid cyclogenesis in 
the Pacific around 150°W. No warm current is present in this region. However. 
Roebber maintains that the maximum is still due to increased baroclinity associated 
with outbreaks of Arctic air through the Bering Straits and across warmer water. 

Since low-level baroclinity is important to rapid cyclogenesis over the open 
ocean, similar conditions may be even more important in coastal regions. The 
proximity of warm water to cold land surfaces produces extremely large horizontal 
temperature gradients in a conditionally unstable environment. Baroclinicallv unstable 
perturbations developing in this favorable environment can be vertically coupled with 
the short waves or jet streaks aloft. 
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D. SUPERPOSITION 



Sanders (1986) showed that rapid intensification takes place as a 500 mb vorticity 
center is superposed to the west of the incipient surface cyclone. The relative position 
between the developing surface low and the vorticity center changes during the life- 
cycle of the cyclone. In composites of rapid deepeners, Sanders showed that the 
average distance between the 500 mb vorticity maximum and the surface low pressure 
center decreases and the vorticity center moves from northwest to southwest of the 
surface cyclone and then to the south as the low develops. 

Farrell (1985) used an initial value approach to show that jet-level and surface 
features exhibit significant growth when properly superposed. An important 
requirement for rapid growth is that there exist a "Fourier overlap" between the upper 
and lower-level disturbances. That is, the upper-level forcing must have significant 
energy in the horizontal scale that matches the surface low. If the superposition 
between relatively large-amplitude features of similar horizontal scale is favorable, 
robust growth can occur even when the waves are neutral or stable to exponential 
growth. 

In summary, rapid cyclogenesis is a phenomenon that depends on interactions 
between processes at all levels of the troposphere, and on several spatial scales. A 
single dominant process during rapid intensification has not been identified, and it may 
vary in different conditions. Additionally, a unique factor that differentiates rapid 
cyclogenesis from "normal" cyclogenesis has not been established. The ERICA 
hypothesis (rapid cyclogenesis requires the proper superposition of PVA aloft and an 
existing surface cyclone in a favorable environment) represents current scientific 
understanding of the phenomenon. This study investigates the nature of the "proper 
superposition" factor. The importance of superposition is first discussed in terms of a 
frictionless, adiabatic model that can be solved analytically. 
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III. TWO-LEVEL MODEL ANALYTICAL STUDY 



A. THE TWO-LEVEL MODEL 

To examine the most basic interactions between low-level and upper-level 
perturbations, an analytical study is conducted using a simple two-level model. The 
quasi-geostrophic equations in pressure coordinates are combined to form a finite 
difference potential vorticity equation following Haltiner and Williams (1980). 

This simplified linear model provides a method to investigate different 
combinations of upper-level and lower-level perturbations in different mean conditions. 
The static stability and the vertical shear in the mean flow (which is directly 
proportional to the meridional temperature gradient) can be changed to study the 
impact on the development of waveform disturbances. The wavelength, amplitude and 
phase relationship of the initial upper and lower perturbations can also be varied. 

By using the two-level model, the most basic baroclinic instability growth 
mechanism can be investigated. Increased growth rates are expected for decreases in 
static stability and increases in mean shear, which can be related to increases in jet 
strength. By inputting finite amplitude waves of different magnitudes and phases at the 
two levels, the relative growth rates can be studied. The effects of an intense jet streak 
or shortwave aloft moving over a weak surface system can be estimated and compared 
to the growth which results when the same upper-level feature is superposed over a 
relatively vigorous surface cyclone. Finally, the effects of different "alongstreanT 
superpositions can be estimated by altering the phase relationship between the upper 
and lower perturbations. 

The phase relationship between upper and lower-level perturbations is described 
intuitively by Hoskins et al. (1985), in which the authors argue logically that the 
induced perturbation velocities at the two levels will reinforce each other, which keeps 
the perturbations in step and enhances growth. They also show that the perturbations 
are locked on to each other with a phase difference that depends on the spatial scale of 
the disturbances. If the phase is changed slightly, the amount of reinforcement also 
changes, which alters the relative phase speeds and forces the perturbations back to 
their equilibrium phase relationship. The two-level model is used to investigate this 
behavior for a relatively small synoptic-scale disturbance in a low stability 
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environment. The interaction between jet-Icvel perturbations and lower-troposphere 
cyclonic centers is an important consideration in the study of rapid cyclogenesis. 
Without a mechanism to vertically connect the features, mutual intensification is not 
possible. 
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Figure 3.1 The vertical structure of the two-level model. 

The vertical * structure of the model is, shown in Fig. 3.1.. The atmosphere is 
divided into four layers, with the wind specified at the intermediate levels (1 and 3) and • 
vertical motion specified at the middle level. For this study, p $ = 1000 mb and 
Ap = 400 mb. 

Assuming waveform solutions for the streamfunction at the upper and lower 
levels, and for vertical velocity in the middle troposphere, the frequency equation for 
the system can be determined. The streamfunctions at levels l and 3 are: 
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and the phase speed is: 
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where p = wave number, 

X = 2fQ 2 '(Ap) 2 <T-) is a stability parameter, 
fQ = Coriolis parameter, 

Ap = pressure interval (py - pQ = 2Ap), 
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a 2 = static stability at level 2, and 
Uj 3 = basic flow at levels 1 and 3. 



The shortwave cutoff is defined by the relation k = ^i c 2 . For ji greater than n c 
(shorter wavelength) the phase speed will be real and the wave will be neutral. The 
system will be unstable when ji 4 < k 2 , provided (Uj-U-j) 2 is sufficiently large. In this 
case, the phase speed can be written as c = c ± ic- where 

* M 



_ U) + _ <3(2/i +A) 

r “ 2 2 /x 2 (/z 2 +A) 



(3.3) 



and 



. „ WXU-UQ*- xVj 1/2 

1 2 ix 2 (jj. 2 +X) 

% 

Then (3.1) can be rewritten: 
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To simplify the mathematics involved in solving the two-level model, two 
assumptions are made. First, the equations are solved on an f-plane (P = 0). Second, 
the mean wind is set equal to zero by making Uj +- U-j = 0. This condition is 
satisfied by making the wind at level 1 westerly and the wind at level 3 easterly with the 
same speed. Then c f = 0 and 
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The streamfunction field at level 1 can be related to that at level 3: 
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Assuming sinusoidal initial conditions for the streamfunctions, the coefficients 
a + ’ and a" are: 
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(3.9) 
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where \}/j 0 and q/^ 0 are the amplitudes of the streamfunctions at time zero. 

Inserting equations (3.6) through (3.9) into (3.5), the streamfunctions at levels 1 
and 3 can be determined for any combination of wave number, basic flow and static 
stability. Under the p = 0 assumption, any wave longer than the shortwave cutofT 
will be unstable, provided L'j is not equal to zero. 

By combining the vorticity equations at levels 1 and 3 and averaging over x, an 
eddy kinetic energy (EKE) equation can be formed: 
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The first law of thermodynamics can be used to form an equation for the change 
in eddy available potential energy (EAPE): 



d(EAPE) 

dt 
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The right side term of (3.10) is found with the opposite sign on the right side of 
(3.11). This term is the conversion of EAPE to EKE. Since the (\}/j - vj/^) term is 
proportional to the temperature perturbation, this term is the two-level model variation 
of the co'T' baroclinic instability conversion of EAPE to EKE. 
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Thus, for a direct circulation (warm air rising and cold air sinking), EAPE is converted 
to EKE. 

A second term contributes to the change in EAPE. (Uj-U^) is the shear in the 
mean flow, and is proportional to the meridional temperature gradient. The x- 
derivative of the eddy streamfunction can be expressed as the meridional wind, v'. 
Thus, this term is proportional to -<3T,'<3y(v'T') and represents the conversion from 
zonal available potential energy (ZAPE) to EAPE. 
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For a linear model such as this, the mean flow is assumed to be constant. 
Therefore, the (Uj-U-j) shear remains constant and provides an unlimited source of 
ZAPE. If the perturbation is unstable, this ZAPE is continually converted to EAPE, 
and then is available for conversion to EKE. Thus unstable perturbations grow 
indefinitely in this model. 



B. BAROCLINICALLY UNSTABLE CASES 

The validity of the model is first tested against the theories of baroclinic 
instability in a standard atmosphere. The static stability was chosen to correspond to a 
constant lapse rate of approximately 6.25 K km' 1 . For this study, the growth rate, 
exp (fic-t), can be written 



Growth rate = exp 
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For a shear of 20 ms* 1 between 800 mb and 400 mb, the maximum growth rate is 2.38 
day' 1 at a wavelength of 4000 km (Fig. 3.2). The shortwave cutoff for this static 
stability -is at 2591 km. 

The growth of perturbations in a standard atmosphere is investigated in 
Appendix A. The model faithfully reproduces expected results. When a 4000 km 
perturbation is inserted, it grows more than fivefold in the first 48 hours and shifts so 
that a westward tilt with height is established. For a 2600 km perturbation (slightly 
longer than the shortwave cutoff) a similar pattern develops, but the growth rate is 
noticeably smaller. The importance of the shear in the mean flow is also verified in 
Appendix A. As expected from (3.14), the growth rate is directly related to the shear. 
When the shear is doubled with a 4000 km perturbation, fivefold growth occurs in less 
than 24 hours. 

As maritime cyclogenesis is generally a small synoptic-scale phenomenon which 
occurs in a low stability environment, the unstable cases of interest are run using a 
static stability lower than that of the standard atmosphere (equivalent to a lapse rate of 
8 K km' 1 ). As the static stability decreases (A increases), the growth rate increases and 
the wavelength of maximum growth shifts to shorter waves, as can be seen by 
comparing the low stability growth rate to the standard atmosphere growth rate in 
Fig. 3.2. The maximum growth rate occurs near 2800 km, and this wavelength is 
chosen for the perturbations. 
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GROWTH RATE 




Figure 3.2. Growth rate for a standard atmosphere (lapse rate 6.25 K km and for^a 
low stability environment (lapse rate 8 K km' 1 ) when the mean flow shear is 20 m s . 

The amplitudes of the initial streamfunction perturbations at levels 1 and 3 are 
varied in Cases 1 through 3 to simulate the superpositions of strong upper features 
over weak or strong surface cyclones. Similarly, weak upper features in Cases 4 and 5 
are superposed over strong and weak lower-level waves. In all of these cases, the shear 
in the mean flow is 20 ms' 1 between levels 1 and 3. with westerly wands aloft and 
easterlies near the surface. 

Cases 1 through 3 investigate the differences in development of low-level 
perturbations which have differing initial amplitudes. The initial amplitude of the 800 
mb streamfunction in Case 1 is 1 nrs' 1 . The corresponding upper-level perturbation is 
an order of magnitude larger (10 nrs' 1 ). The initial perturbations are in phase (there is 
no tilt with height). The ensuing growth of both perturbations is shown in Fig. 3.3. 
The striking development in this experiment is the extremely rapid growth at the lower 
level. While the 400 mb perturbation triples in size in 36 hours, the amplitude of the 
low-level streamfunction increases from 1 nrs' 1 to 32.5 nrs' 1 . As will be seen, this 
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Perturbation Amplitude 

Case 1 




Figure 3 . 3 . Amplitudes of Case 1 400 mb (\ j / j , solid) and S00 mb (>1/3, dashed) 

• perturbations through 36 h. 

rapid growth to an amplitude nearly as large as that at the other level is a 
characteristic of every case in which one of the initial amplitudes is much smaller than 
the other. 

For Cases 2 and 3 the level 1 perturbation is again 10 nrs* 1 . Case 2 is initiated 
without a perturbation at the lower level. Case 3 investigates an initial lower 
perturbation which is as large as the upper feature. As in Case 1, the initial 
perturbations are in phase. The growth of the 800 mb perturbations for Cases 2 and 3 
is compared to that for Case 1 in Fig. 3.4a. The corresponding phase relationships are 
shown in Fig. 3.4b. A positive phase difference indicates a westward tilt with height. 

Although Case 2 is initiated without an 800 mb perturbation, a low-level feature 
still develops and the pattern is not appreciably different from case 1. After only 9 
hours the upper feature lags the lower by approximately 70° in both cases. The low- 
level amplitude in Case 2 actually exceeds that from Case 1 after 6 hours. These 
results indicate that the presence of a weak closed circulation at the surface does not 
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CASES 1,2, and 3 

Lower Perturbation Amplitude 





Figure 3.4. (a) Amplitudes of 800 mb perturbations, and (b) phase differences between 
400 mb and 800 mb perturbations through 36 h for Cases 1 (solid), 2 (dotted), 3a 
(dashed), and 3b with an initial westward tilt of 90° (heavy dashed). 
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appreciably enhance development when a strong in-phase jet-level feature is 
superposed. Indeed, a weak surface low in the wrong position (Case 1) will actually be 
detrimental to development compared to a case without an initial low-level 
perturbation (Case 2). 

In Case 3 the initial low-level perturbation is as large as the upper feature. The 

amplitudes at both levels after 36 hours are somewhat larger than for the cases with 

little or no initial low-level disturbance, although the overall intensification at the lower 

level is actually smaller than in cases 1 and 2. As neither level dominates in Case 3, the 

wave grows steadily at both levels. There is no period of extremely rapid growth as the 

smaller wave grows to the same order of magnitude as that of the other level, as there 

was in Case 1. After only 9 hours in that case, the lower perturbation amplitude was 

more than five times the initial value. Additionally, the phase difference between the 

two levels after 9 hours is nearly a quarter-cycle. On the other hand, when the lower 

perturbation is as large as at the upper level (Case 3) the increase in the low-level 

amplitude early in the experiment is not nearly as dramatic. The growth in 9 hours is 

less than half that experienced in Case 1. Also, the phase separation between the 
) 

upper and lower waves early in the period .is much smaller. 

It might be expected that the Case 3 growth, which starts with large 
perturbations at both levels, should far exceed the growth in cases where the low-level 
perturbations are small. However, Case 3 growth is significantly smaller than that 
expected (.3.4 day' 1 -) from the growth rate equation (3.14). The Case 3 perturbation, 
actually takes 23 hours to double in size. The energetics equations (3.12) and (3.13) 
help explain this apparent discrepancy. 

The conversions from ZAPE to EAPE and EAPE to EKE are proportional to the 
correlations between the eddy quantities v', co' and T'. The sizes of these individual 
terms are important, but the energy conversions also depend on the spatial 
relationships. Improper phasing can result in an indirect circulation (cold air rising as 
warm air sinks) which means EKE is converted to EAPE and the perturbation is 
actually decaying. 

Synoptic experience indicates that a perturbation which is vertically stacked is 
not favorable to further development. This "barotropic" state is most often associated 
with the mature and decaying stages of synoptic-scale cyclones. A more favorable 
pattern for baroclinic growth is a perturbation which tilts westward with height, 
"against" the mean flow shear. It can be seen in Fig. 3.4b that this favorable tilt is 
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quickly established in Cases 1 and 2. However, it takes 15 hours for a 45° phase 
difference to develop in Case 3. This lack of favorable superposition explains the 
failure of the Case 3 perturbation to grow as rapidly as expected. 

Cases 1 and 3 are repeated with an initial westward tilt of a quarter-wavelength. 
In both cases growth is faster, but the difference is especially noticeable in Case 3 
(heavy dashed line in Fig. 3.4). The growth rate in this case becomes that predicted by 
(3.14). The perturbation doubles in amplitude in only 13 hours because the 
superposition is favorable for baroclinic growth from the beginning of the experiment. 

In each of these cases the phase difference clearly approaches the same value 
(approximately 67°). This westward tilt with height is the equilibrium phase difference 
described by Hoskins et al. (1985). For larger (smaller) phase differences, the total 
perturbation velocity causes the disturbances to move toward (away from) each other 
until the equilibrium tilt is established. 

For Case 4, the size of the initial upper and lower perturbations are reversed, 
with the 10 irrs'^ amplitude at 800 mb and the 1 irrs'^ amplitude at 400 mb. The 
growtlrat the lower level (Fig. 3.5) is very similar to Case 1, even though the initial 
perturbation at the upper level is small. If the large perturbation is initially at the 
lower level, the spectacular growth is at the upper level, and the amplitude quickly' 
approaches that of the 800 mb wave. As in Case 1, the tilt of the wave after 36 hours 
is toward the west with height, with a phase of about 67°. This is to be expected in 
the two.-level model, since friction is neglected and the static stability is assumed- 
constant with height. The equations for the streamfunctions at levels 1 and 3 are 
symmetric, so if the initial conditions are simply inverted, with no other changes to the 
mean flow or stability, the result must be growth rates which are also reversed. 

The development of the perturbations in Cases 1 and 2 can be thought of as 
Petterssen Type B events, (Petterssen and Smebye, 1971) in w r hich an upper-level 
positive vorticity maximum is superposed over a region which is favorable for 
baroclinic instability, but w'hich has little or no initial vorticity at the surface. The 
Case 4 development does not fit this pattern. A large perturbation at the surface under 
an upper flow which is nearly zonal causes extremely rapid growth at the upper level. 
This case is more like a Petterssen Type A disturbance, in w r hich self-development 
causes upper-level amplification as the surface disturbance grows and creates 
differential temperature advection. Resulting thickness changes cause an upper-level 
ridge; trough system which is situated such that positive vorticity advection aloft will 
cause further intensification of the surface low. 
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Figure 3.5. Amplitudes of 400 mb and 800 mb perturbations through 36 h for Case 4. 

Vj / 1 (solid) and (dotted); and for Case 5, V 3 (dashed). 

If both the upper and lower-level initial waves are small, growth is very small 
compared to cases in -which either or both of the features is large. In case 5, the 
amplitudes of both initial waves are set at 1 mV 1 . Although the phase shift between 
the growing waves is similar to the rapid growth cases, the growth is much slower than 
in cases 1 through 4. This is not surprising, since the horizontal and vertical velocity 
and temperature perturbations are all quite small. Thus, the conversions from ZAPE 
to EAPE and EAPE to EKE remain small and subsequent growth is limited. 

It can be concluded from these five cases that growth through baroclinic 
instability will result when any perturbation is inserted in a mean flow which contains 
sufficient vertical shear (horizontal temperature gradient). The magnitude of the 
growth depends on the size of the initial perturbations. If the eddy is large at either 
400 mb or 800 mb (or at both levels), subsequent growth is large. These Finite-size 
perturbations can be thought of as linear model representations of observed synoptic 
features commonly associated with rapid cyclogenesis. Actual cases of cyclogenesis are 
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rarely, if ever, attributed to the growth of small-amplitude perturbations in flow which 
is initially undisturbed. Rather, cyclogenesis is typically associated with the coupling of 
large-amplitude jet-level perturbations (shortwaves and jet streaks) with low-level 
cyclonic circulations in a favorable environment. 

The two-level model indicates that the vertical location of the initial large 
perturbation is not critical to the resulting growth. Similar intensification will result if 
the initial large wave is at either 400 mb or 800 mb, which implies that a strong surface 
wave under steady flow aloft can induce rapid cyclogenesis as effectively as a 
propagating jet streak or upper-level shomvave. The addition of friction and vertically- 
varying static stability to the model may alter this conclusion, but the two-level model 
alone does not favor upper-level processes as more effective in inducing cyclogenesis. 

C. PHASE DIFFERENCES IN THE INITIAL PERTURBATIONS 

Two other sets of experiments are made with the Case 1 and 3 initial conditions 
to further study the effect of phase differences between the initial perturbations. Using 
the Case 1 initial conditions, the small low-level initial feature was displaced relative to 
the large upperdevel feature and the development was studied for a full range of phase 
differences (Fig. 3.6). The initial conditions can be written 

V|/j(x,0) = 10 cos px , (3.15) 



and 



vy^x.O) = cos p(x + 5) . (3.16) 

where 6 is the fraction of a wavelength phase difference between the two perturbations. 

When Cases 1, 3, 4 and 5 from the previous section are integrated with westward 
tilts (90° phase difference) in the initial conditions, the growth rates increase over the 
experiments with the perturbations in phase. A full range of initial tilts are tested to 
determine the most favorable and least favorable superpositions. In all cases, the 
phase differences converge to 67° (Fig. 3.6b). The only appreciable differences are in 
the amplitudes of the perturbations, where the fastest growing case (5 = 115°) exceeds 
the least favorable case (5 = -65°) by 20 percent after 36 hours. 
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Figure 3.6. fa) Amplitudes of 800 mb perturbations and (b) phase differences in case 1 
through 36 h with initial phase differences of 6= -135° (dashed), 6= -65° (dotted). 
6= 45° (dot-dashed), and 6= 115° (solid). 
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A closer look at the growth during the early part of the period is useful. In every 
case the phase of the low-level perturbation quickly shifts so that a westward tilt is 
established. Beyond 18 hours the phase difference is approximately 67° for every case. 
As expected from baroclinic instability theory, maximum growth occurs when the low- 
level wave is initially placed in a westward tilt position. When the 800 mb perturbation 
leads the 400 mb wave by 90° to 135°, growth in the first several hours is significantly 
greater than when the initial perturbations are unfavorably phased. However, after the 
first few hours, the phase relationship does not appear to be very important for a large 
perturbation superposed over a small one. Regardless of the initial positioning, the tilt 
becomes favorable for baroclinic growth within six hours. The energy concentrated at 
400 mb dominates the growth, quickly establishing a lower-level perturbation at a 
favorable position. All perturbations then intensify rapidly, with the growth in the 
least favorable case only four hours behind that in the most favorable case. 

It is hypothesized that the initial phase relationship between the upper and lower 
features will be much more important when the perturbations are approximately the 
same size (Case 3). Neither wave dominates, and it takes longer for a shift to a 
favorable westward tilt with height to be established. In this - case the level 3 
perturbation can be written 

Vj/^x.O) = cos)i(x + S). (3.17) 

Again the initial low-level wave is displaced relative to the upper-level perturbation. If 
the tilt is "downshear" the waves are actually damped during the early period (Fig. 3.7). 
If the initial tilt is westward with height, the waves at both levels grow from the outset. 
After 9 hours, the amplitude of the waves in the most favorable case has grown nearly 
65 percent over the initial amplitude, while the amplitude in the least favorable case 
has decreased 37 percent. 

In contrast to the small amplitude case in Fig. 3.6b, the phase difference changes 
slowly during the early development when large amplitude perturbations are present at 
both levels (Fig. 3.7b). As the upper perturbation induces an 800 mb disturbance 
downstream, the lower perturbation induces a 400 mb wave upstream. The actual 
pattern is a summation of the two, and the disturbances gradually move to the 
equilibrium position. In the least favorable case (5 = -65°) the phase difference 
remains negative through the first 44 hours. For other initial phase differences, the tilt 
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(a) Amplitudes of 800 mb perturbations and (b) phase differences in Case 
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Figure 3.7. 

through 36 h with initial phase differences of 5= -135° (dashed). 6= -65° (heavy 
dashed). 6= -45 (dotted). 8 = 45° (dot-dashed), and 8 = 115° (solid). 
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becomes westward with height somewhat sooner. When the orientation becomes 
favorable for growth, the wave begins to grow at both levels. The amplitudes after 36 
hours are markedly different, which demonstrates the importance of the initial 
superpositioning when both features are comparatively large. If the initial tilt is 
favorable (6 = 90° to 8 = 135°), the 36 hour growth is nearly sevenfold, but the 
perturbation continually decays to less than 15 percent of the initial value for an initial 
tilt of -65°. 

The most favorable initial phase difference is 115° for both Case 1 and Case 3. 
This superposition represents a balance between the conversion of EAPE to EKE 
(greatest at 180°) and the conversion of ZAPE to EAPE (maximized at 90°). As the 
phase difference shifts towards 67°, this case experiences the greatest early growth. 
The 67° phase difference is where the change in EAPE (conversion from ZAPE minus 
loss to EKE) is the most positive. 

It is concluded that the importance of the initial perturbation superpositioning 
depends on the relative size of the features. If the energy at the jet-level greatly 
exceeds that near the surface, the lower-level wave will quickly shift/redevelop so that a 
westward tilt with height is established. If both the jet-level and the surface 
circulations are strong, the superpositioning is critical to further development. In this 
case, the upper wave cannot quickly overpower the low-level feature and establish a 
favorable westward tilt. If the initial phasing is unfavorable, the upper and lower 
waves interfere and the amplitudes at both levels, decrease. Eventually a westward tilt 
is established and growth begins through baroclinic instability. By this time, the 
magnitude of the unfavorably phased perturbation is significantly smaller than that of 
the wave which had a favorable tilt at 'the initial time. Therefore, the efficiency of the 
energy transfers from ZAPE to EAPE and EAPE to EKE will be higher for the 
favorable tilt case. After 36 hours, the amplitude of the 8 = 115° case perturbation is 
nearly 50 times the amplitude of the 8'= -65° wave. Comparing Fig. 3.7a to Fig. 3.6a, 
when the Case 3 large amplitude perturbations are favorably placed, the lower-level 
growth is significantly greater than that for the most favorably phased Case 1 example 
with a large perturbation aloft over a small perturbation at the surface. 

D. NEUTRAL WAVE DEVELOPMENT 

The static stability in the two-level model is constant with height and in time. It 
is recognized that the value used in Cases 1 through 5, which corresponds to a lapse 
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rate of 8 K km' 1 , may be unrealistically low for an average through the entire 
troposphere. Although low stability values are often associated with oceanic 
cyclogenesis (Sanders and Gyakum, 1980; Bosart, 1981; Sardie and Warner, 1983; and 
many others), these values are usually confined to the lower layers of the atmosphere. 
Nevertheless, the wavelengths of many explosively deepening cyclones 'are much 
shorter than the 4000 km predicted for an average mid-latitude atmosphere. 

Recent papers by Farrell (1982, 1984, and 1985) and Shepherd (1985) discuss the 
short-term transient growth of waves which are baroclinically stable. The development 
of neutral waves is also investigated using the two-level model, and the results are 
presented in Appendix A. For the standard atmosphere static stability, the shortwave 
cutoff is 2591 km. In this f-plane study, any perturbation longer than the cutoff is 
unstable. Several perturbation wavelengths shorter than the critical value are tested. 
The neutral waves do not maintain a favorable vertical connection between the 
perturbations. The phase continually changes throughout the experiment, and the 
amplitudes of the neutral waves grow (decay) when the phase difference is positive 
(negative). When the wave tilts "against" the shear, the conversion from ZAPE to 
EAPE to EKE continues in a manner similar to baroclinic instability. However, the 
phase will eventually shift to an eastward tilt. As soon as the phase difference changes 
sign, the amplitude of the wave begins to decay. The closer the wavelength is to the 
critical value, the longer the period before decay begins and the larger the maximum 
amplitude of the perturbation. 

Different initial phases are also tested on a neutral case (wavelength 2000 km). - 
As in the unstable cases, introduction of initial phase differences when one 
perturbation is much larger than the other does not greatly affect subsequent 
development. Low-level perturbations in the most favorable case (5 = 180°) reach a 
maximum amplitude only 2 nr s' 1 greater than in the least favorable case. 

In the unstable cases, initial phase differences were critical when the upper and 
lower features were the same size. In an analogous neutral case the initial tilt of the 
perturbation is again crucial to the development of the wave (Fig. 3.8). The least 
favorable tilt for a neutral wave is when the initial perturbation is vertically stacked. 
As can be seen in Fig. 3.8b, the neutral wave continually shifts toward a more eastward 
tilt. If the perturbations are initially in phase, the phase difference immediately 
becomes negative and the wave converts EKE to EAPE and continues to decrease in 
amplitude until the phase difference passes ISO degrees. At that time, the next wave in 
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Figure 3.8. (a) Amplitudes of 800 mb perturbations and (b) phase differences through 
96 h for a neutral case (L = 2000 km) with initial phase differences of 6= 0° (dashed). 
0= 90° (dotted), and 6= 180° (solid). 
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the upper-level wave train is in position to support growth and the perturbations grow 
back to their initial size. 

The most favorable superposition for the neutral case occurs when the initial 
perturbations are 180 degrees out of phase. In this case, the early development results 
in wave growth of nearly 100 percent before the phase difference becomes negative and 
the amplitude again decreases. The maximum amplitude of the oscillating perturbation 
in the most favorable case is double that produced when initial superposition is least 
favorable. 

E. CONCLUSIONS 

A two-level linear analytical model is used to investigate basic relationships 
between upper-level and lower-level perturbations. While this type of linear model is 
too simple for operational use, it has great utility for the study of the most basic 
relationships between waves at an upper level (approximating the jet level) and 
perturbations present above the planetary' boundary layer. 

The energy transformations through which the potential energy of the basic state 
(evident in the shear between levels 1 and 3) is fed to the kinetic energy of the 
perturbations can be equated to the baroclinic instability transfers first identified by 
Charney (1947) and Eady (1949). Since the model is linear, the -eddies do not interact 
with the mean flow, nor is the available potential energy of the mean flow decreased 
during integration. An infinite source of energy is therefore available for exponential 
growth of perturbations at both levels. 

As long as perturbations remain relatively small, the linearity assumption is 
reasonable and variations in growth rates among the tested cases are indicative of what 
might be found in actual events. Since the model is linear, the size of the initial 
features can be scaled as desired to ensure that they remain comparatively small 
throughout each integration. The major interest here is not the actual amplitude of the 
perturbations, but the relative size of the growth rates at the two levels and for 
different cases. 

Many of the conclusions reached using the two-level model are as expected from 
theories derived from more complicated models and observational studies. For 
example, energy transformations and perturbation growth rates are significantly larger 
when static stability is decreased, and the growth rate is directly proportional to the 
shear between levels 1 and 3. Baroclinic instability causes growth whenever the wave 
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tilts westward with height (for mean flows which become more westerly with height). 
That is, the growing waves "oppose" the shear. Study of shorter neutral waves reveals 
that significant growth can take place before the tilt first becomes eastward with 
height. 

Several two-level model conclusions are important to the study of explosive 
cyclogenesis, and will be further examined in following sections: 

• When a relatively large perturbation is superposed over a small one (or even a 
lower level which is. without perturbation), the lower level perturbation 
experiences a period of extremely rapid growth. 

• The phase relationship between the initial perturbations is critical if the 
perturbations are about the same size at both levels. If the superposition is 
such that the wave tilts eastward with height, the early development will be 
damped at both levels until a favorable tilt is established. When the 
perturbations are favorably superposed, growth through baroclinic instability 
can begin at the outset. Superposition of the upper and lower features is also 
important for the neutral case if both perturbations are relatively large. The 
phasing is not very important if the energy is primarily at one level. This level 
will dominate the other and a favorable westward tilt will quickly be established. 

• For a given static stability and mean flow, maximum growth occurs when a 
relatively large upper-level perturbation is favorably phased with respect to a 
well-developed surface feature. 
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IV. NORAPS CHANNEL MODEL EXPERIMENTS 



A. THE MODEL 

The conclusions reached in the two-level model study of Chapter III are further 
tested using a numerical model that provides a more realistic representation of the 
atmosphere. The model chosen is the Navy Operational Regional Atmospheric 
Prediction System (NORAPS), which is described in detail by Hodur (1982, 1984). 
NORAPS is a primitive equation limited area model in <7 coordinates. Momentum and 
mass variables are placed on a type C staggered grid (Arakawa and Lamb, 1977). The 
model uses split-explicit time integration, fourth-order advection, and fourth-order 
diffusion. The grid is globally relocatable with variable horizontal resolution. The 
number and location of the <7 surfaces can be chosen by the user. One-way interactive 
boundary conditions are provided from the corresponding N'OGAPS (Navy 
Operational Global Atmospheric Prediction System) run. The full model includes a 
Deardorff(1972) planetary boundary’ layer scheme, dry convective adjustment, large- 
scale precipitation, and a modified Kuo-tvpe cumulus parameterization (Kuo. 1965). 

The experiments in this chapter use a channel version - of NORAPS to make the 
numerical experiments analogous to the analytical study of Chapter III. For these 
experiments, the grid spacing is set at 100 km in the x and v directions. The channel is 
6400 km wide with free-slip, impermeable boundaries on the northern and southern 
sides and cyclic boundary’ conditions in the east-west direction. The channel length 
depends on the wavelength of the longest wave in each experiment. There are 12 sigma 
levels in the vertical between the surface and 50 mb. Frictionless and adiabatic options 
are chosen, and p is set to zero. 

Because the experiments in Chapter III were made with a linear model. 
NORAPS is run in a quasi-linear mode in this chapter. The effects on the mean flow 
due to nonlinear interactions with the perturbations are removed by restoring the mean 
flow after each step of the integration. That is. the zonal average winds, temperatures 
and heights at each level are replaced with the initial zonal mean values of u.v.T and <p 
during each time step. This method allows development and movement of all 
perturbations without changing (reducing) the initial sources of zonal available 
potential energy (ZAPE) and zonal kinetic energy (ZKE). Wave-wave interactions do 
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occur, but are assumed to be negligible, especially for perturbations longer than global 
wave number 20 (Hodur, 1984). 

A preliminary growth-rate experiment is conducted to determine the most 
unstable wavelength for conditions that parallel those of the two-level model cases. 
The mean wind is a zonal jet stream that is centered in the middle of the channel (Fig. 
4.1). The width of the stream is exaggerated in the central 2000 km- wide swath to 
reduce any possible effects of shear vorticity and to require a constant north-south 
temperature gradient. This imitates the conditions in the two-level model, which did 
not include variations in the y-direction. The vertical shear in the central region is 
constant with height from U = 0 at the surface to U = 45 m s' 1 at 200 mb. The winds 
above 200 mb are constant with height. This wind profile makes the shear between 
400 mb (U 400 = 27 m s' 1 ) and 800 mb (U 800 = 7 m s' 1 ) equal to the 20 m s* 1 used in 
the two-level cases. The mean temperature lapse rate (8 K km' 1 ) is equivalent to the 
low-stability lapse rate of Chapter III. 

For the growth-rate experiment, initial sinusoidal perturbations are inserted in a 
103 x 65 grid for model wave numbers 1 through 10 (wavelengths 10,000 km down to 
1000 km). The perturbations consist of two half-sine waves in cr-coordinates so there is 
a double maximum in amplitude in the vertical (Fig. 4.2). In' this experiment, the 
initial amplitudes of the upper and lower perturbations are the same. The upper 
maximum is centered at 400 mb (as in the two-level model) and extends from 187 mb 
to 613 mb. The lower maximum is centered at 800 mb (again paralleling the two-level 
study) and extends from -588 mb to the surface. These perturbations have a sine- 
squared profile in the cross-channel direction, with zero amplitude at the northern and 
southern walls and maximum amplitude at mid-channel. 

After the perturbations are inserted, the model is initialized and integrated for 
84 h. Growth rates (Fig. 4.3) are determined for each wave number by analyzing the 
change in the v-component with time. This curve is reassuringly similar to the low 
stability growth rate curve for the two-level model (Fig. 3.2). The wavelength of 
maximum growth is near 2800 km (as in the two-level model), and the maximum 
growth rate is 3.28 day' 1 , which compares to 3.38 day* 1 in the analytical study. 
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Figure 4.1. Cross-section of the mean zonal wind (U) field for NORAPS channel 
model Cases 1 through 5. 
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Figure 4.2. Vertical and cross-channel shape of initial upper and lower perturbations. 



50 



GROWTH RATE 




Figure 4.3. Growth rates observed in the NORAPS channel model for lapse rate 
S K km' 1 and the wind profile of Figure 4.1. 

3. BASIC NORAPS EXPERIMENTS 

1. Experiments Without Initial Phase Differences 

Based on the growth rate curve, a 2S00 km perturbation is chosen for the 
remaining NORAPS channel model experiments. The x-dimension of the channel is 
reduced to one wavelength to save integration time. Thus, each experiment includes 
only wavelengths equal to or shorter than the most unstable perturbation (model wave 
number one). Experiments analogous to all Chapter 3 cases are made with the 
NORAPS channel model (Table 1). As in the analytical study, an order of magnitude 
difference in initial amplitude is used to differentiate between strong and weak 
perturbations. In the two-level model, the perturbations were inserted as small changes 
to the streamfunction field. In this chapter, the initial waves are introduced as v- 
component winds. The case numbers are the same as in Chapter 3 for ease in 
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comparison, except that Case 4B in this chapter is an additional case. The magnitudes 
of the upper and lower-level v-prime values in Table 1 (0, 1 and 10 cm s' 1 ) are the 
maximum amplitudes of the upper and lower perturbations. As they did in the 
analytical study, the maxima occur in the center of the channel at 400 mb and 800 mb. 
Since the mean flow is purely zonal, the v-field at any time can be used to establish the 
position and amplitude of the developing perturbations. 



Case 


TABLE 1 

NORAPS Channel Model Cases 
Upper v' (cm s' 1 ) Loner v' (cm s' 1 ) 


1 


10 


l 


2 


10 


0 


3 


10 


10 


4A 


1 


10 


4B 


0 


10 


5 


1 


1 



Cases 1 through 3 explore the response of weak, zero and strong lower 
perturbations under strong upper-level forcing. Cases 4A and'4B test weak and zero 
upper waves with strong low perturbations, while weak waves at each level are included 
in Case 5. As in the first set of two-level cases, these perturbations are inserted at the 
same zonal position. That is, there is initially no vertical phase difference between the 
perturbations. 

In all cases, the perturbations develop similar profiles of perturbation 
amplitude versus height (Fig. 4.4). The peak values (10 cm s' 1 ) at the initial time are 
at 400 and 800 mb. Around 600 mb, the initial perturbations overlap and the 
amplitude is a sum of contributions from the upper and lower features. The only 
growth during the first 12 h of the integration is at the surface as the perturbation is 
re-distributed among the levels. By 36 h. the basic profile of the growing wave is 
evident. The maximum amplitude is at the surface, as expected for a frictionless model. 
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Figure 4.4. Case 3 perturbation amplitude versus height. Initial conditions (bold). 

T = 12 h (dotted), t = 24 h (dashed), z = 36 h (dot-dashed), t = h (thin solid). 



and a secondary maximum is found near 380 mb. A similar profile develops in every 
case, even when the initial perturbation is absent at the lower (Case 2) or upper level 
(Case 4B). If the lower perturbation maximum amplitude had been originally placed at 
the surface rather than 800 mb, the 24 h adjustment period would have been much 
shorter, and the amplitudes of the resulting low-level perturbations would have been 
slightly larger. However, the relative growth of the waves in the different cases would 
have been unchanged. 

The phase of the perturbation at each level is determined by locating the 
maximum v-component and expressing the distance from a fixed reference point as a 
phase angle, where 360° equals one wavelength (2800 km). As expected, a westward 
tilt with height is quickly established and is maintained for the duration of the 
experiment (Fig. 4.5). The plots of phase versus height can be quite irregular during 
the first 18 to 24 h of the experiments. Again, had the initial low-level amplitude 
maximum been placed at the surface, these oscillations would have been reduced. For 
comparison with the two-level study, the phase differences between 400 mb and 800 mb 
are calculated by fitting a least-squares line to the phases at all levels below 350 mb. 

The first conclusion from the two-level study (that lower perturbations will 
experience a period of rapid growth when a vigorous upper feature is superposed) is 
verified by these NORAPS integrations. In all frictionless NORAPS cases, the 
amplitude of the lower- level (specifically, at the surface) perturbation exceeds that at 
the upper level after 24 h of integration, even if the initial conditions contain no lower- 
level perturbation. The development through 48 h of the Case 3 perturbations at both 
levels is summarized in Fig. 4.6. As in the two-level study, a westward tilt with height 
is quickly established and the phase difference eventually reaches an equilibrium value 
that is maintained for the remainder of the experiment. This preference for a certain 
phase difference is described intuitively by Hoskins et al. (1985) and was summarized 
in Chapter III. Unlike the two-level study, the 400-800 mb phase difference exhibits a 
damped oscillation as the perturbations approach equilibrium separation. As the upper 
and lower perturbations adjust toward positions where they provide maximum positive 
feedback, they overshoot these positions and are then forced back until eventually a 
balance is reached. In the two-level model, the equilibrium phase separation between 
400 and 800 mb was approximately 67°. In all NORAPS cases, the eventual 
equilibrium phase difference is 55°. If the phase difference between 400 mb and the 
surface (where the low-level perturbation is largest) is calculated, the most favorable 
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Figure 4.5. Vertical tilt of Case 3 perturbations through 4S h. 
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Figure 4.6. (a) Amplitudes of upper (between 200 and 500 mb) and lower (below 700 
mb) perturbations and (b) phase differences between 400 mb and 800 mb through 48 h 
for Case 3. Positive (negative) phase differences indicate westward (eastward) tilt with 
height. 
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separation is nearly 73°. Therefore, the equilibrium separation that is eventually 
established is nearly the same in both models. 

The amplitudes of the initial perturbations in the NORAPS study are 
relatively large compared to that of the streamfunction perturbations in the two-level 
study. Since linear theory assumes infinitesimal perturbations, the validity of finite- 
amplitude waves in this "quasi-linear" model may be suspect. If the quasi-linear 
approach is valid, the amplitudes at any time should be proportional to the amplitude 
of the initial perturbations. Using Case 3 as a basis for comparison, additional 
integrations are made in which the initial amplitudes are multiplied by 10 (v' = 
1 m s' 1 ) and 50 (v' = 5 m s* 1 ). The upper and lower perturbation amplitudes for the 
largest case are shown in Fig. 4.7. Except for the scale of the y-axis, this case almost 
exactly parallels the small amplitude Case 3 (Fig. 4.6a). Therefore, the quasi-linear 
approach is satisfactory for initial perturbations of this size. 

Cases 1, 2 and 3 combine a vigorous upper perturbation with small, zero and 
large lower features, respectively (Fig. 4.8). Only minor differences exist between these 
cases and the analogous two-level model cases (Fig. 3.4). In the NORAPS simulations. 

y 

the Case 1 amplitude remains larger than Case 2 throughout the integration. Whereas 
the Case 1 and 2 amplitudes at 36 hours are only about 50 percent of that in Case 3, 
the Case 1 and 2 perturbations in the two-level study grew to about 90 percent of the 
amplitude of Case 3. In the tw’o-level model, the Case 3 phase difference shifted to the 
equilibrium separation in about 30 h. Using NORAPS, the Case 3 tilt becomes nearly 
the same as the Case 1 tilt after about 12 h. 

In summary, similar relationships exist in both models. When both 
perturbations are relatively strong (Case 3), initially unfavorable phase differences in 
both models are slower to shift to the equilibrium separation compared to Cases 1 and 
2. In the two-level study, the Case 3 phase difference was slow' to reach the optimal 
westward tilt, which resulted in growth that was slow r er than expected. During the first 
16 h of the N’ORAPS Case 3 integration, the perturbation tilt is eastward or slightly 
westward and the amplitude of the lower perturbation remains nearly constant during 
this period of adjustment. In either model, rapid growth occurs when a favorable 
westward tilt is established. One difference between the two experiments is that the 
zero initial phase difference superposition is not as unfavorable for growth in NORAPS 
as was indicated by the two-level model. 
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Figure 4.7. Amplitudes of upper and lower perturbations when Case 3 initial 
perturbation amplitudes are multiplied by 50. 
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Figure 4.S. (a) Amplitudes of low-level perturbations and (b) phase differences through 
4S h for Cases 1 (dashed), 2 (dotted), and 3 (solid). 
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The lower perturbation amplitudes and phase differences for the remaining 
cases are shown in Fig. 4.9. In Case 4A, an initially strong lower perturbation is 
superposed below a weak upper perturbation. In Case 4B, the upper level is without an 
initial perturbation. The developments in Cases 4A and 4B are similar to Cases 1 and 
2, with the small perturbation case (4A) growing only slightly faster than the no- 
perturbation case (4B). Case 5 initially contained small perturbations at both levels. 
Without large values in the eddy wind, mass and temperature fields, the energy 
conversions in Case 5 are small compared to those of Cases 1 through 4 and the net 
growth of the perturbation is very small. 

2. Effect of Initial Zonal Phase Differences 

An important conclusion from the tw'o-level study is that the initial zonal 
positions of upper and lower-level perturbations are critical when both perturbations 
are relatively strong. For example, initial perturbations that were vertically stacked 
(zero initial phase difference) proved to be unfavorable for growth. In the previous 
section, it was indicated that the in-phase position did not -seem to be as detrimental to 
Case 3 growth in NORAPS. Some differences between the tw r o models that might 
contribute to differences in the NORAPS study include;. 

• numerical versus analytical solutions; 

• "quasi-linear" versus linear approximations; 

• vertical resolution that is greatly increased and vertically-varying static stability 
in NORAPS; and 

• mass and temperature fields that are fixed in the two-level model are free to 
adjust to each other in NORAPS. 

To test the effects of changes in the initial phase difference with NORAPS, the 
upper perturbation (10 cm s' 1 in Cases 1 and 3, 1 cm s* 1 in Case 4A) is placed at the 
same location for each experiment, and the lower perturbation is moved east and west 
along the centerline in the channel until the entire range of possible zonal phase 
relationships is tested. Based on the two-level model, the effect of the initial phase of 
the superposition is minimal when one perturbation is much larger than the other 
(Cases 1 and 4A). In the NORAPS cases, the phase differences also quickly shift to a 
westward tilt in all cases and converge to nearly the same value by 9 h. The differences 
in growth of the lower wave (Fig. 4.10) for Case 1 (upper v’ = 10 cm s' 1 , low r er v' = 

1 cm s' 1 ) are largest between initial phase differences of 6 = 55° (most favorable) and 
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Figure 4.9. As in Figure 4.S, except for Cases 4A (solid), 4B (dashed), and 5 (dotted). 
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Figure 4.10. (a) Amplitudes of low-level perturbations and (b) phase differences 
through^ h for Case 1 with initial phase differences of 6 = -135° (heavy dashed) 
6 - -4o (dashed). 5 = 0° (dot-dashed), 6 = 55° (solid), and 5 = 180°'(dotted). 
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5 = -135° (least favorable). The amplitude is reduced by about 25 percent in the least 
favorable case. While the equilibrium phase separations that are eventually reached 
are nearly the same in the NORAPS and two-level models, the initial superpositions 
that result in the greatest (and smallest) growth are quite different. In the two-level 
study (Fig. 3.6), the most (least) favorable initial superposition was 6 = 115° (-65°). 
The slightly larger response to the initial tilt in the NORAPS model may be due to the 
tendency for even more rapid low-level growth in NORAPS than in the two-level 
model. As the lower perturbation grows (even in the 5 = -135° case), it induces 
upper-level development to the west, W'hich interferes with the existing large upper-level 
perturbation and slows the establishment of the most favorable westward tilt. In Case 
4A (small upper perturbation superposed on a large low-level feature) the effect of 
changing the initial phase difference is even less important. The amplitudes at both 
levels in the least favorable case are only about 12 percent less than the most favorable 
case (not shown). 

Where strong features exist at both levels (Case 3) in the two-level model, the 
phase difference is very slow to shift to a favorable westward tilt if it is initially 
unfavorable (Fig. 3.7)/ In contrast, the phase differences in the NORAPS Case 3 
display extremely rapid shifts (Fig. 4.11b) and the upper and lower perturbations 
continually overshoot the equilibrium value. Thus, the phase differences are' still slow 
to converge, which has the same effect on the amplitudes as in the two-level study. 
That is, cases that start near the equilibrium value (6 = 55°) grow quite rapidly, since 
the large perturbations are positioned in such a way as to allow maximum energy 
conversion from zonal available potential energy (ZAPE) to eddy available potential 
energy (EAPE) and from EAPE to eddy kinetic energy (EKE). On the other hand, the 
phase difference in the least favorable case (5 = -135°) oscillates wildly for nearly 36 h 
before finally converging on a 55° westward tilt. Consequently, the perturbations are 
in unfavorable superpositions as often as they are favorably positioned, and the 
amplitude of the 6 = -135° wave at 48 h is about equal to the initial amplitude. 

Energy budget terms for selected Case 3 runs are shown in Table 2. The 
energy budget equations are those of Robertson and Smith (19S3) as modified for the 
closed domain by Hodur (1984). The dominant budget terms are the baroclinic 
instability processes: CA— conversion of zonal available potential energy (ZAPE) to 
eddy available potential energy (EAPE); and CE— conversion of EAPE to eddy kinetic 
energy (EKE). The next most important terms are the barotropic term (CK), which 
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Figure 4.11. As in Figure 4.10, except for Case 3. 
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TABLE 2 



Energy Conversions for Case 3 in units of 10' 3 W m' 2 . See tex,t 
for explanation of terms. 



Tilt 


Time ( h ) 


CA 


0 ° 


12 


0.15 




24 


0.65 




36 


2.56 




48 


9.14 


55 ° 


12 


0.15 




24 


0 . S2 




36 


■ 3.09 




48 


10.66 


O 

o 

CO 


12 


0.03 




24 


0.25 




36 


0 . 67 . 




48 


2.10 


135 ° 


12 


0.03 




24 


0.08 




36 


0.11 




48 


0.45 


- 45 ° 


12 


0.11 




24 


0.37 




36 


1.47 




48 


5 . 4S 



CE 


CZ 


CK 


o.os 


0.01 


0.07 


0.42 


0.03 


0.02 


1.66 


- 0.07 


0.10 


5.34 


- 0.02 


0.22 


0.13 


0.0 


0.01 


0.57 


0.02 


0.02 


1.94 


- 0.14 


0.08 


-6.17 


- 0.07 


0.21 


0.05 


0.0 


'o.o 


0.15 


0.0 


0.0 


0.39 


- 0.05 


- 0.02 


1.20 


- 0.05 


0.0 


0.01 


0.01 


- 0.01 


0.01 


0.01 


0.0 


0.09 


0.03 


- 0.01 


0.31 


0.0 


0.0 


0.03 


0.02 


0.0 


0.20 


0.02 


0.02 


0.99 


0.0 1 


0.07 


3.25 


0.01 


0.16 
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measures the conversion of zonal kinetic energy (ZKE) to EKE, and the conversion 
(CZ) of ZAPE to ZKE by vertical redistribution of mass. However, the magnitudes of 
CK and CZ are only about 5 percent of the leading terms. The baroclinic terms are 
positive at all times (instability) and generally increase exponentially through the 
integration. However, the magnitudes of the leading terms vary widely for different 
initial superpositions. For example, CA and CE in the least favorable case 
(6 = -135°) are only about 5 percent of the most favorable case (5 = 55°). 

Thus, both the linear model and NORAPS results demonstrate that the initial 
phase difference is crucial when upper-level and lower-level perturbations are both 
strong. The most favorable initial tilt differs between the two models and the phase 
difference changes slowly in the two-level model and rapidly in NORAPS. While 
shifting from a 115° phase difference to the equilibrium position, the two-level model 
perturbations grow rapidly as a result of consistently large baroclinic energy transfers. 
In NORAPS, the same initial superposition results in oscillations around the 
equilibrium phase difference that are larger than the oscillations that occur when the 
initial phase difference equals the equilibrium value. While the most favorable initial 
tilts are different in the two models, both agree in general with synoptic' Experience and 
theoretical studies (Eadv, 1949; etc.) that a westward tilt with height of about a 
quarter-wavelength provides optimal conditions for intensification. 

3. Summary of Zonal Tilt Results 

The use of a quasi-linear, frictionless numerical model with a constant (low) 

static stability throughout the troposphere confirms the three most important 

conclusions reached in the two-level analvtical studv. 

✓ ^ 

• When a large perturbation is present at the jet level and. there is little or no 
energy at lower levels, a low-level perturbation will be induced and will undergo 
a period of extremely rapid growth. When this low-level perturbation is small 
relative to the upper-level perturbation, the initial phase of the superposition is 
almost irrelevant. That is, the growth of a new perturbation near the most 
favorable position relative to the large upper-level perturbation dominates over 
the contribution from the initial low-level perturbation. 

• For a given mean flow and static stability, a more rapid growth occurs when 
the perturbation at the lower level is also relatively large and the initial 
superposition leads to a westward tilt of 55° with height. 
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• The initial phase difference of the superposition is critical when both upper and 
lower perturbations are nearly the same strength, with maximum (minimum) 
growth when the phase difference is 55° (-135°). 

Thus far, only along-stream (zonal tilt) superposition has been investigated. 
The NORAPS model will be used to examine the cross-track superposition in the next 
section. 

C. EFFECT OF CROSS-CHANNEL SUPERPOSITION 
1. Procedure 

The superposition questions investigated in the previous sections were 
restricted to along-channel (zonal) tilts. In every case, the upper and lower 
perturbations were centered under a broad jet stream in the middle of the channel. 
Real-world disturbances are not confined to a position under the jet stream. It is 
common for cyclogenesis to occur south of the jet stream and for the intensifying 
system ‘to cross under the jet at some point in its development. The average position 
of rapidly deepening surface cyclones is in the cold air to the north of the associated 
500 mb forcing (Sanders and Gyakum, 1980; Rogers and Bosart, 1986) but this 
relationship is not as well defined as the westward tilt with height. Other studies have 
found rapidly deepening cyclones located south of (or under) the jet stream (Jannuzzi, 
1985; Sinclair and Elsberry, 1986). The added dimension of a cross-track tilt will be 
investigated in this section. 

To study the effect of y- variations in the perturbation locations, it is necessary 
to use a more realistic mean flow than that from the previous sections (Fig. 4.1). The 
jet stream shown in Fig. 4.12 has the same vertical structure in mid-channel as 
previously, in that the speed increases linearly with height from 0 ms' 1 at the surface 
to 45 m s' 1 at 200 mb and remains constant above this level. To satisfy thermal wind 
constraints, the horizontal temperature gradients at each level in mid-channel are also 
equivalent to those in previous sections. For example, the surface temperature 
gradient is approximately 1.1 K, 100 km. The mean winds north and south of the jet 
decrease symmetrically to the boundaries so that horizontal shear effects are 
introduced. 

The perturbations shown in Fig. 4.2 are modified slightly by making the y- 
dependencv a function of sin 4 y instead of sin 2 y, which concentrates perturbation 
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Figure 4.12. As in Figure 4.1, except for the "basic" case cross-track tilt experiment. 
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energy nearer the central latitude of the wave. The upper-level perturbation is still 
centered at the mid-point of the channel, but the lower-level feature is placed north or 
south of mid-channel as well as east or west of the upper wave. The Case 3 amplitudes 
(10 cm s' 1 at both levels) are used in this section, since growth is most sensitive to 
superposition when both perturbations are relatively strong. The perturbations are 
placed in a zonal superposition that is favorable for growth (55° westward tilt from 800 
mb to 400 mb) and the lower wave is positioned between 800 km north and 800 km 
south of the jet. 

Many factors may influence the most favorable cross-track position for the 
surface low pressure center. Not all of them can be included in this frictionless, 
adiabatic model that does not include along-stream variations in the mean flow (jet 
streaks). Two factors that can be investigated are the relative contributions of the 
surface baroclinity and shear vorticity effects to growth of the low-level perturbation. 
After the first set of experiments using the basic jet (Fig. 4.12), tests are conducted 
using a mean wind with enhanced boundary layer shear (strong surface baroclinity) and 
using a mean wind profile with enhanced mid-tropospheric shear. A fourth experiment 
is conducted with a mean flow jet- that slopes from south to north with height rather 
chan being vertically stacked as in Fig. 4.12. The four- cases are designated "basic." 
"baroclinic," "shear" and "tilt" for ease in reference. In all cases, the mean lapse rate 
(8 K km’ 1 ) is' equivalent to that used in previous experiments. 

2. Basic Jet Case 

Perturbations centered at 400 mb and 800 mb are inserted in the mean flow 
shown in Fig. 4.12 with the upper disturbance in the center of the channel (north- 
south) and the 800 mb wave placed 55 0 to the east at a variety of positions (-800. -600. 
-400, -200, 0, 200, 400. 600, 800 km) to the north and south. Further integrations are 
then made to isolate the most favorable meridional tilt. After 24 h. the largest 
amplitude perturbations are found at the surface in every integration. Based on the 
relative sizes of these surface disturbances at 48 h (Fig. 4.13), the most favorable 
meridional superposition is slightly north (60 km) of the center of the channel for this 
mean flow. The maximum growth rate is 1.62 d" 1 . The differences in 48 h amplitude 
among the different north-south positions are relatively small since the zonal phase 
difference is favorable for growth in every case. The growth in the least favorable case 
(800 km south) is only 20 percent less than that for the 60 km north case. Experiments 
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Figure 4.13. Amplitudes of surface perturbations after 4S h for cases with low-level 
perturbations initially displaced from mid-channel. Basic case (solid), baroclinic case 
(dot-dashed), shear case (dashed), tilt case (dotted). 
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with the initial low-level perturbation more than 1000 km from the most favorable 
meridional position indicate that the growth rate decreases even more rapidly with 
increasing separation as the perturbation is removed from the strongest tropospheric 
baroclinic zone. 

With the addition of y-tilt in the problem, the most favorable east-west tilt 
may be other than 55°. Several combinations of zonal and meridional tilts in the 
vicinity of 55°, 60 km north are tested, and it is found that this combination results in 
maximum growth for this particular mean flow. 

The perturbations are of course free to move meridionally and zonally relative 
to the initial locations. Just as the zonal phase differences gradually converged on the 
equilibrium value of 55°, the lower perturbations in these integrations move toward the 
most favorable north-south position as they grow (Fig. 4.14). After 48 h. the average 
positions of the maximum and minimum surface pressure centers for all the basic jet 
cases are within 150 km of the most favorable position at 60 km north of mid-channel. 
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Figure 4.14. Latitudinal positions of surface disturbances after 48 h for the basic case. 



The leading energy budget terms for several of the experiments in this case are 
presented in Table 3. The magnitudes of all energy transfers for initial positions 
greater than 400 km north or 600 km south of the center of the channel are smaller 
than for the initial positions nearer mid-channel. The baroclinic instability terms CA 
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TABLE 3 

Conversions for the Basic Jet Case 


in units of 10' 4 


W m’ 2 . 


Tilt 


Time (h) 


CA 


CE 


CZ 


CK 


400 km N 


12 


0.87 


0.43 


0.02 


-0.02 




24 


3.47 


2.36 


0.03 


0.08 




36 


10.23 


6.34 


-0.36 


0.45 




48 


29.90 


16.19 


-0.24 


1.48 


200 km N 


12 


0.89 


0.49 


0.03 


-0.01 




24 


3.4S 


2.35 


0.05 


0.08 




36 


10.42 


6.42 


-0.33 


0.46 




48 


30.31 


16.56 


-0.25 


1.46 


60 km N 


12 


0.88 


0.54- 


0.03 


0.0 




24 


3.48 


231 


0.06 


0.08 




. 36 


10.48 . 


6.42 


-0.32 


0.44 




48 ' 


30.25 


16.57 • 


-0.23 


1.39 


200 km S 


12 


0.84 


0.59 


0.03 


0.0 




24 


3.37 


°2.15- 


0.10 


0.07 




36 


10.19 


6.14 


-0.31 


0.39 




48 


28.96 


15.97 


-0.21 


1.23 


400 km S 


12 


0.79 


0.62 


0.03 


0.01 




24 


3.22 


2.01 


0.13 


0.05 




36 


9.73 


5.79 


-0.28 


0.33 




48 


27.16 


15.05 


-0.19 


1.05 


600 km S 


12 


0.74 


0.62 


0.03 


0.01 




24 


3.01 


1.82 


0.16 


0.03 




36 


9.03 


5.30 


- 0.27 


0.25 




48 


24.75 


13.71 


-0.16 


0.84 




and CE are by far the largest energy conversions, with magnitudes more than 20 times 
as large as the barotropic term (CK) and the conversion of ZAPE to ZKE (CZ). At 
the earliest time period (12 h), the largest conversion from ZAPE to EAPE (CA) is 
slightly north of the most favorable position, while the largest CE (conversion from 
EAPE to EKE) actually occurs for perturbations 500 km south of mid-channel. 
However, the perturbation at 60 km north grows fastest due to a combination of all 
the energy transfers. As this perturbation continues to grow larger than those in the 
other experiments, it naturally has the largest magnitude energy, conversions for all 
times greater than 12 h. 

3. Enhanced Surface Baroclinity Case 

One factor that may be expected to be extremely important in rapid 
cyclogenesis is surface baroclinity. Enhanced boundary layer baroclinity has been 
suggested as an ingredient in a number of cases of rapidly intensifying cyclones 
(Sanders and Gyakum, 1980; Bosart and Lin. 1984). A mean flow that features a 
surface temperature gradient of 1.95 K/100 km (approximately 75 percent greater than 
that in the basic case) is similarly tested for the most favorable cross-channel tilt. The 
maximum velocity of the jet is maintained at 45 m s' 1 , but the vertical wind shear 
between 850 mb and 200 mb is reduced. Thus, there is a more rapid increase in mean 
wind velocity between the surface and 850 mb than in the basic jet case (Fig. 4.15). 
Thermal wind considerations require that this vertical shear be balanced by a large 
temperature gradient in the lowest levels. 

The meridional position of the fastest growing wave in this "baroclinic" case 
shifts to 390 km south of mid-channel (Fig. 4.13). As in the basic case, all 
perturbations grow since the zonal tilt is westward with height. However, the amount 
of growth by 48 h decreases more rapidly with distance from the most favorable 
position at 390 km south. For example, the growth in the 800 km north case is 50 
percent less than the maximum growth. The 23 percent increase in the maximum 
growth rate (2.0 d' 1 ) over the basic jet case is attributable to the increased ZAPE in the 
baroclinic jet (Table 4). In every case, the perturbations converge toward a position 
about 500 km south of mid-channel by 48 h (Fig. 4.16). As in the basic case, the 
southern part of a disturbance placed to the north of the most favorable latitude will 
grow fastest, which shifts the perturbation slightly to the south. This process continues 
until the wave reaches the latitude of maximum growth. 
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Figure 4.15. As in Figure 4.1, except for the "baroclinic" case cross-track tilt 
experiment. 
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TABLE 4 



Total ZAPE (10 5 J m‘‘) and maximum growth rates (d' 1 ) during the first 
48 h for the cross-channel cases. 



Case 


ZAPE 


Maximum Growth Rate 


Basic 


45.6 


1.62 


Baroclinic 


57.3 


2.00 


Shear 


.51.2 


1.64 


Tilt 


47.5 


1.51 



Barocilnic Case 

Y-Positlon of surface perturbations 




The sensitivity to small zonal and meridional displacements about this 55° 390 
km south position in the "baroclinic" case was tested. The maximum growth rate 
occurs for a superposition in which the low-level perturbation is initially 45° (350 km) 
east of and 390 km south of the upper-level wave. However, this slightly smaller zonal 
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phase difference results in only a 1 percent increase in growth rate over the 55°/390 km 
south case. 

Changes in the mean flow may slightly alter the most unstable wavelength. 
Particularly, increased low-level baroclinity may decrease the wavelength of maximum 
growth. The experiments in this study were all conducted with 2800 km perturbations. 
The results of experiments with a shorter wavelength should show only minor 
differences, and the major conclusions should remain the same as in the 2S00 km case. 

The southward shift in the most favorable position in the baroclinic case may 
be due to an enhanced frontogenesis mechanism when the surface temperature gradient 
is initially strong. The surface frontal zone forms to the south of a zonal jet stream 
(Reed and Danielson,. 1959). Stronger frontogenesis south of the jet results in greater 
baroclinic energy conversion in that region. Since vertical motions and temperature 
perturbations are especially enhanced in frontal regions, the conversion from EAPE to 
EKE (CE), w r hich is proportional to -to'T', should be especially important in the region 
of the front. The leading energy terms for a number of integrations in this case are 
showrn in Table 5. The energy terms for cases with initial perturbations farther north 
are all smaller than those in the table. Early in the integrations (12 h) when the 
perturbations are still near their initial positions, the maximum CA is in the 390 km 
south case and the maximum CE is in the 800 km south case. The importance of the 
coT' contribution can best be illustrated by the ratio of CE to CA at 12 h. The farther 
south the perturbation is placed, the more important CE becomes. The ratio increases 
° from 0.32 for the 800 km north case to 0.71 for the 800 km south case. By 48 h, all the 
perturbations are south of the jet and the CE to CA ratio is 0.47 for every case, as 
expected for a linear system, after the most unstable mode has become dominant. 

Thus, the most favorable position (390 km south) represents a balance 
between the case with the most rapid early growth and the case with the largest energy 
growth later in the integration. At 12 h, the 100 km south case has the largest 
amplitude. After 36 h, the 600 km south case has the largest conversion to EKE, w’hich 
indicates that this perturbation w r ill eventually dominate. 

4. Enhanced Upper-Tropospheric Shear Case 

Another important factor in the meridional location of developing cyclones is 
the relative vorticity present in the mean flow. For this experiment, the relative 
vorticity of the mean wind is due only to horizontal shear, since the mean wind field 
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TABLE 5 

Energy Conversions for the Baroclinic Case in units of 10' 4 W m' 2 . 



Tilt 


Time (h) 


CA 


CE 


CZ 


CK 


200 km N 


12 


1.44 


0.72 


-0.03 


-0.03 


• 


24 


5.34 


2.98 


-0.01 


0.05 




36 


18.60 


8.98 


-0.84 


0.21 




48 


60.76 


28.57 


-1.08 


0.21 


Center 


12 


1.51 


0.84 


-0.05 


-0.02 




24 


5.74 


3.11 


0.01 


0.06 




36 


20.32 


9.66 


-0.85 


0.19 




48 


67.04 


31.41 


-1.07 


-0.01 


200 km S 


12 


1.56 


0.94 


-0.05 


0.00 




24 


6.10 


3.20 


0.03 


' 0.06 




36 


21.83 


10.28 


-0.86 


0.14 




48 


72.82 


34.12 


-1.03 


-0.33 


390 km S 


12 


1.58 


1.01 


-0.05 


0.01 




24 


6.31 


3.21 


' 0.05 


0.04 




36 


22.76 


10.62 


-0.84 


0.06 




48 


76.42 


35.80 


-0.94 


_ -0.69 


600 km S 


12 


1.57 


1.06 


-0.05 


0.01 




24 


6.37 


3.18 


0.06 


0.03 




36 


23.16 


10.75 


-0.81 


• -0.04 




48 


78.15 


36.62 


-0.86 


-1.15 


800 km S 


12 


1.53 


1.08 


-0.06 


0.02 




24 


6.26 


3.06 


0.08 


0.01 




36 


22.84 


10.53 


-0.78 


-0.15 




48 


77.22 


36.20 


-0.74 


-1.55 



77 



has no curvature. The effect of the enhanced shear vorticity is tested by increasing the 
mean wind above 500 mb by a factor up to 1.33 in the central third of the channel. 
The new maximum wind is 60 m s' 1 at mid-channel from 200 mb to the top of the 
model (Fig. 4.17). Since the vertical wind shear in the lower troposphere is unchanged, 
the surface temperature gradient and low-level baroclinity are the same as in the basic 
jet case. This increased shear may result in barotropic instability at upper levels south 
of the jet. Although the magnitude of the barotropic energy conversion term (CK) is 
larger than in previous cases, it is still only 10 percent of the baroclinic terms, CA and 
CE (Table 6). As in the baroclinic case, the changes in the mean flow may slightly 
alter the most unstable wavelength. 

As in the previous cases, an upper-level perturbation is inserted at mid- 
channel in this "shear" case and lower perturbations are tested in a variety of 
meridional positions relative to the jet stream. All integrations result in reasonable 
growth, regardless of the initial meridional position of the lower perturbation (Fig. 
4.13). However, the greatest growth occurs in the case in which the perturbation is 
initially 200 km north of the lower jet. That is, the increased middle and upper-level 
shear results in a northward shift of the most favorable postion. Further sensitivity 
tests indicate maximum growth for this mean flow is achieved when the superposition 
includes a 50° westward tilt and a 200 km north displacement. 

The positions of the surface perturbations after 48 h demonstrate pronounced 
northward shifts in these "shear" cases (Fig. 4.18). In every, case, the surface 
perturbations converge toward a position approximately 350 km north of the jet 
stream. As in the previous experiments, perturbations initially placed near this position 
have the largest energy conversions at 48 h (Table 6). Early in the integrations (12 h) 
the largest conversion from ZAPE to EAPE (CA) coincides with the most favorable 
position at 200 km north of the jet for this experiment. This case is also similar to the 
baroclinic case in that the location with the largest amplitude at 48 h is between the 
location with the greatest early growth (largest perturbations at 12 h are for the mid- 
channel case) and the location with the maximum energy conversion at later times. 

Increasing the jet strength also results in a small (less than 5 percent) increase 
in the maximum lower perturbation growth rate compared to the basic case. This 
increase in the intensity of the surface low is smaller than might have been expected, 
since the ZAPE for this case is 12 percent greater than in the basic case (Table 4). 
This increased ZAPE is introduced in the upper troposphere and the relatively small 
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Figure 4.17. As in Figure 4.1, except for the "shear" case cross-track tilt experiment. 
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TABLE 6 



Energy Conversions for the Shear Case in units of 10* 4 W m' 2 . 



Tilt 


Time (h) 


CA 


600 km N 


12 


0.90 




24 


4.72 




36 


14.73 




48 


45.62 


400 km N 


12 


0.93 




24 


4.73 




36 


14.96 




48 


45.90 


200 km N 


12 


0.94 




24 


4.70 




36 


14.94 




48 


45.25 


Center 


12 


0.91 




24 


4.54 




36 


14.49 




48 


43.26 


200 km S 


12 


0.86 




24 


4.30 




36 


13.74 




48 


40.22 


400 km S 


12 


0.79 




24 


3.95 




36 


12.60 




48 


36.11 



CE 


CZ 


CK 


0.50 


0.09 


0.00 


2.92 


0.12 


0.20 


8.76 


- 0.76 


1.11 


25.39 


i 

o 

o 


4.54 


0.56 


0.10 


0.01 


2.93 


0.12 


0.20 


8.89 


- 0.77 


1.13 


25.68 


-0.50 


4.58 


0.56 


0.11 • 


0.02 


2.90 


0.13 


0.20 


8.84 


- 0.76 


1.10 


25.41 


-0.08 


4.49 


0.64 


0.11 


0.03 


2.78 


0.15 


0.18 


8.53 


-0.75 


1.03 


24.31 


-0.10 


4.26 


0.65 


0.12 


0.05 


2.61 


0.16 


0.16 


S.04 


-0.71 


0.92 


22.67 


-0.10 


3.88 


0.63 


0.12 


0.06 


2.39 


0.17 


0.13 


7.35 


-0.68 


0.80 


20.37 


-0.06 


3.44 
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Figure 4.18. As in Figure 4.14, except for the shear case. 



vertical shear below 500 mb is not conducive to transport of the increased energy 
toward the surface. It is also surprising that the magnitudes of the upper perturbations 
(as measured by the v-component of the wind) do not increase relative to the basic 
case. This may be a reflection of the reduced vertical energy flux, or it may be a 
feature of the quasi-linear approach, in which the mean zonal wind is restored to the 
original, value after each time step. Although the. upper level v' values are essentially 
the same, the northward flux of zonal momentum is much larger for the enhanced jet. 
In summary, increasing the jet velocity without increasing lower tropospheric 
baroclinity tends to shift the most favorable position for surface low development to 
the north of mid-channel, but does not result in significantly larger perturbation 
growth. 

5. Tilted Jet Case 

The final experiment in this section is conducted with a mean flow jet that 
slopes northward with height (Fig. 4.19). Although the winds above 200 mb are the 
same as the basic case, the meridional location of the maximum wind speed tilts from 
the center of the channel at 200 mb to 1000 km south of mid-channel at the surface. 
This 1,100 slope is more realistic than the symmetric cross-sections used in the 
previous cases (Reed, 1955; Reed and Danielson. 1959). The maximum surface 
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Figure 4.19. As in Figure 4.1, except for the "tilt" case cross-track tilt experiment. 
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temperature gradient is the same as in the basic case, but it is centered 380 km south of 
the center of the channel. 

As in the other cases, the westward tilt makes all positions favorable for 
growth. The meridional position of the initial perturbation that results in the largest 
amplitude at 48 h is 230 km south of the jet (Fig. 4.13). This position is under the 
sloping jet core between the maximum surface temperature gradient and the positive 
shear vorticity. Based on the "baroclinic" case results, the southward shift of the 
baroclinic zone should force a similar shift in the meridional position of the most 
favorable growth case. However, the upper-tropospheric shear is nearly the same as 
the basic case, which favors surface low intensification just north of mid-channel. The 
relatively large separation between the tilts favored by the shear and baroclinitv factors 
results in growth- rates in this "tilt” case that are the smallest (maximum 1.51 d* 1 ) of 
any in this section, even though the ZAPE in the mean How is slightly larger than in 
the basic case (Table 4). The most favorable case experiences net growth that is 41 
percent greater than for the least favorable case (800 km north). The largest 
conversions to EKE early in the integration is for the case with the perturbation at 600 
km south (Table 7). After 12 h, the largest baroclinic energy conversions are for the 
perturbation that starts just north of mid-channel. In the tilted case, the perturbations 
have a decreased tendency to move toward the most favorable latitude compared to 
previous experiments (Fig. 4.20). The conflicting influences of the low-level baroclinitv 
.to the south and positive relative vorticity aloft to the north reduce the lateral 
movement of the perturbations. After 48 h, the surface pressure centers are still 
distributed over a range of 785 km, whereas the separation is only 270 km in the basic 
case (Fig. 4.14). The superposition that results in the largest amplitude perturbation at 
48 h is a compromise between the latitude of the greatest early growth (600 km south) 
and the latitude of the wave that will eventually dominate (50 km north). 

6. Summary' of Cross-Channel Tilt Results 

Several conclusions can be drawn from these experiments in which different 
meridional tilts are tested within a simple model (frictionless, adiabatic, and quasi- 
linear). 

• When upper-level and lower-level perturbations of approximately the same 
amplitude are superposed, the most important factor determining subsequent 
growth is the along-stream (east-west) tilt between them. Changes in the 
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TABLE 7 



Energy Conversions for the Tilt Case in units of 10‘ 4 W m" 2 . 



Tilt 


Time (h) 


CA 


CE 


CZ 


CK 


400 km N 


12 


0.70 


0.33 


0.07 


-0.03 




24 


2.36 


1.77 


0.15 


0.04 




36 


• 6.22 


3.93 


-0.12 


. 0.31 




48 


17.00 


9.31 


0.03 


0.97 


200 km N 


12 


0.71 


0.38 


0.08 


-0.03 




24 


2.39 


1.75 


0.16 


0.04 




36 


6.40 


3.97 


-0.11 


0.31 




48 


17.34 


9.57 


0.04 


0.94 


Center 


12 


0.71 


0.43 


0.08 


-0.02 




24 


2.41 


1.70 


0.18 


0.04 




36 


6.44 


3.91 


' -0.11 


0.29 




48 


17.28 


9.53 


0.06 


0.87 


230 km S 


12 


0.68 


0.47 


0.07 


-0.01 




24 


2.39 


1.59 


0.21 


0.03 




36 


6.35 


3.78 


-0.12 


0.25 




48 


16.66 


9.22 


0.09 


0.73 


400 km S 


12 


0.65 


0.49 


0.07 


-0.01 




24 


2.33 


1.50 


0.23 


0.02 




36 


6.18 


3.60 


-0.12 


0.21 




48 


15.97 


8.82 


0.11 


0.60 


600 km S 


12 


0.62 


0.51 


0.07 


0.00 




24 


2.27 


1.40 


0.25 


0.01 




36 


5.95 


3.36 


-0.14 


0.16 




48 


15.01 


8.29 


0.13 


0.46 
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Figure 4.20. As in Figure 4.14, except for the tilt case. 
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cross-flow (north-south) position by less than 1000 km only slightly modify the 
growth rates established by the zonal superposition. 

• Lower perturbations that are laterally displaced from the most- favorable 
latitude shift toward this position as they grow. The adjustment in zonal tilt to 
the most favorable value (all cases converged to the equilibrium tilt of 55° by 
48 h) appears to be more rapid than the convergence in the meridional tilt. 

• The most favorable cross-track superposition depends on the characteristics of 
the mean flow. Increased low-level baroclinity enhances the conversion of 
EAPE to EKE south of the jet, and moves the most favorable position of the 
lower wave toward the south. Enhanced horizontal shear in the middle and 
upper-troposphere (increased cyclonic vorticitv north of the jet) moves the most 
favorable position of the lower wave toward the north. 

• The position of maximum 48 h growth in these experiments is a balance 
between the position where early (0-12 h) growth is largest and the position that 
will eventually dominate (greatest energy conversions for times > 36 h). 

• A jet core that slopes from south to north with height results in a region of 
maximum surface temperature gradient that is laterally separated from the zone 
of cyclonic vorticitv. Growth rates in this case are smaller than for a vertically 
stacked jet with the same velocity and surface temperature gradient. 
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D. JET STREAK SUPERPOSITION 



Along-stream maxima in jet velocity (jet streaks) are important factors in upper- 
level forcing of cyclogenesis. Jet streaks provide an additional source of positive 
vorticity advection, and superpose a "checkerboard" convergence/divergence pattern on 
the basic flow. For example, divergence aloft in the left exit and right-entrance 
quadrants can enhance existing upward vertical motion if properly positioned with 
respect to the surface disturbance. Curvature effects may also be added if the jet streak 
is associated with an upper-level shortwave trough. In this section, a jet streak (rather 
than the 400 mb perturbation) is superposed over a finite-amplitude low-level 
perturbation. Growth rates at both levels are investigated for a variety of zonal and 
meridional superpositions. 

.1. Procedure 

A sinusoidal variation in the x-direction with wavelength 2800 km is added to 
the mean flow in the basic jet case (Fig. 4.12). The corresponding v-component is 
derived by specifying that the total flow be non-divergent. A balanced mass field is 
then derived. Although the flow is initially non-divergent and balanced, these 
constraints do not apply during the integration. The resulting jet-level (200 mb) flow 
pattern has a maximum velocity of 52 m s* 1 and a mid-channel minimum velocity of 36 
m s* 1 (Fig. 4.21). As the jet streak provides significant forcing to lower levels, the 400 
mb perturbation that was previously inserted at the beginning of each experiment is 
not necessary. Whereas the width of the perturbation in the y-direction was the same 
at both levels in previous experiments, the y-scale of the jet streak is different from that 
of the 800 mb perturbation. Specifically, v' in the jet streak is maximum 500 km north 
and south of the center of the channel, and is zero along the centerline. This y-scale 
mismatch may slightly reduce the efficiency of the energy transfers to the lower levels. 
As in previous sections, the model is run in a frictionless, adiabatic, quasi-linear mode. 
The quasi-linear approach retains the zonal mean wind throughout the integration, but 
the jet streak, which can be considered an eddy in the zonal mean flow, is free to move 
and develop. 

The first jet streak experiment is conducted without an initial surface 
perturbation. As expected, a surface system is quickly induced by the upper-level 
forcing. This development is used as a baseline value for evaluating further experiments 
that include low-level perturbations at different initial positions with respect to the jet 



86 



y-diatance (km) 



JET STREAK CASES 

200 mb wind (m s' 1 ) 




Figure 4.21. Isotachs (contours at 5 m s" 1 intervals) at 200 mb for the initially 
non-divergent jet streak. 
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streak. In the zonal tilt experiments, the initial 800 mb perturbation is in the center of 
the channel at various zonal superpositions relative to the jet streak. Since the surface 
disturbance induced by the jet streak is significantly larger (about an order of 
magnitude) than the disturbances induced in section C, the amplitude of the initial low- 
level perturbation is increased to 0.5 m s' 1 for these experiments. After the most 
favorable zonal tilt is established, cross-channel superposition with the jet streak is 
investigated. 

2. Jet Streak-Induced Surface Development 

With no initial perturbation at 800 mb, the jet streak quickly forces the 
formation of a surface system below the left-exit quadrant of the streak. The v- 
component amplitudes of the upper and lower-level disturbances and the phase 
difference between the jet streak and the surface trough through 48 h are shown in Fig. 
4.22. As in previous experiments, the v-component of the wind is used as a measure of 
the intensity of the perturbations, and horizontal distances between locations of the 
maximum v-components are used to determine the phase difference. 

This case is analogous to Case 2. in which a large upper-level feature induces , 
perturbation growth in purely zonal- low-level fields. Although the low-level growth 
induced by the jet streak (Fig. 4.22a) is an order of magnitude larger, it follows the 
same general curve as the Case- 2 growth (Fig. 4.8). Upper-level growth is more 
vigorous when a jet streak provides the forcing. Whereas the lower perturbation 
amplitude exceeds the upper after 24 h in Case 2 (not shown), the upper feature 
remains larger through the first 44 h in the jet streak case. The vertical structure of the 
developing disturbance is also slightly different between the two cases. In the jet streak 
case, the phase difference changes rapidly very' early in the experiment (Fig. 4.22b). A 
favorable westward tilt is established after 4 h, and the phase difference converges 
toward an equilibrium separation nearly the same as that in Cases 1 through 5 (Fig. 
4.8b). These phase difference oscillations with a jet streak present are less pronounced 
than in previous experiments. However, the convergence toward the equilibrium 
separation is slower. The separation at 48 h is 82.5°, but is still decreasing toward an 
eventual equilibrium value of 54°. 
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Figure 4.22. (a) Upper and lower perturbation amplitudes and (b) phase differences 
through 4S h for the jet streak case without initial low-level perturbation. Positive 
(negative.) phase differences indicate westward (eastward) tilt with height. 
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3. Jet Streak Zonal Tilt Experiments 

The effect of initial zonal phase differences of a 0.5 m s* 1 lower perturbation 
relative to the jet streak is investigated in this section. Based on the straight jet stream 
experiment and Fig. 4.22, it is expected that a westward tilt with height from the low- 
level perturbation to the streak will be most favorable for development. The phase 
difference is between the center of the jet streak and the 800 mb trough. A full range 
of possible zonal superpositions (0° to 360°)is tested, and the most favorable position 
is identified. 

As in the straight jet stream case, the initial superposition is critical for 
subsequent development. Depending on the initial tilt, the surface perturbation after 
48 h may be larger or smaller than the surface perturbation induced by the jet streak 
alone (Fig. 4.23). The most and least favorable separations (55° and -135°) are the 
same as in the previous case. -Net growth in the -135° case is approximately 40 percent 
less than in the 55° case. The phase differences also follow the pattern established in 
previous sections. After a period of adjustment, a favorable westward tilt is established 
in every case, and the separation between the 800 mb .and 400 mb features gradually 
converges toward the equilibrium value (52°). As in the jet streak case without a 
perturbation, -the phase difference oscillations are smaller and the convergence toward 
the final tilt is slower than in the straight jet case. 

The feedback between the upper and lower levels differs for the jet streak case. 
In previous experiments, the superposition that resulted in the largest surface 
intensification was also most favorable for upper-level growth. However, the greatest 
perturbation growth above 500 mb in this experiment is not in the 55° initial tilt case, 
but when the surface trough and the jet streak are almost vertical (5° westward tilt 
with height ). This nearly "barotropic" stacking places the maximum low-level southerly 
flow approximately a quarter- wavelength ahead of the jet streak, with the left-exit 
quadrant divergence reinforcing the upward vertical motion associated with the low- 
level warm advection. This 5° westward tilt superposition is also reasonably favorable 
for surface intensification, but the resulting lower perturbation amplitude at 48 h is 
smaller than in the case with a 55° initial tilt. 

4. Jet Streak Cross-Channel Tilt Experiments 

The sensitivity of the jet streak case to cross-track position is investigated by 
placing the 0.5 m s' 1 S00 mb perturbation in a favorable zonal tilt (55°) and at a 
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Jet Streak Cases with Initial Phase Differences 
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Figure 4.23. (a) Amplitudes of low-level perturbation and (b) phase differences 
through 4S h for jet streak cases with initial phase differences of 6 = ISO 0 
(dash-dotted), 6 = 55 ° (solid), 6 = 0° (dashed), 6 = -135° (dotted), and no 
perturbation (heavy solid). 
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variety of positions north or south of the jet streak. When the model was run without 
a low-level perturbation, the jet streak induced surface cyclogenesis at a position 
approximately 600 km east and 515 km north of the center of the jet. Presumably, a 
perturbation initially placed near this position would experience the most favorable 
conditions for rapid growth. Indeed, the largest energy conversion terms at 36 and 
48 h are for the 400 and 600 km north cases (Table 8). However, the largest 48 h 
amplitude is found for a perturbation initially placed just south (60 km) of the center of 
the channel (Fig. 4.24). However, the increase in amplitude from the least favorable 
position (800 km north) to the 60 km south position is only 6 percent. 

All of the developing surface centers move to a position approximately 500 km 
north of the jet within the first hour of the integration. Because the favorable position 
is shared by all the perturbations after 1 h, the 48 h amplitude depends directly on the 
amplitude of the perturbations very early in the integration. After the first hour, the 
low-level perturbation in the 60 km south case is largest of all the cross-channel cases 
(approximately 10 percent larger than the 500 km north case). This result again shows 
the importance of the surface baroclinic zone to the very early growth of the 
perturbations. Initial growth is directly related to the low-level temperature gradient, 
which is maximum at mid-channel. The perturbation that starts just south of this 
position is shifted northward during the earliest time steps, benefitting most from the 
strong baroclinity. Maximum growth for the jet stream case is in the 50° westward 
tilt/60 km south superposition, although the amplitude is only marginally larger than 
the 55 O/ 60 km south combination. 

5. Summary of Jet Streak Results 

A non-divergent jet streak is added to the zonal basic flow conditions in the 
previous sections to investigate this form of upper-level cvclogenetic forcing. When the 
model is run without a lower perturbation, the results are analogous to Case 2 (10 
cm s' 1 upper perturbation and zero lower perturbation). The cyclogenesis is much 
more energetic in the jet streak case, but the size of the lower perturbation relative to 
the upper feature, and the phase differences between them are quite similar to Case 2. 
In general, the jet streak experiments support the results from previous sections, with 
only minor modifications. 

• The zonal phase difference in the initial conditions is extremely important to 
subsequent growth. The most favorable initial tilt is 55° between the center of 
the velocity maximum and the initial S00 mb trough. 
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TABLE 8 

Energy Conversions for the Jet Streak Case in units of 10' 2 W m’ 2 . 



Tilt 


Time (h) 


CA 


CE 


CZ 


CK 


800 km N 


12 


1.92 


-0.39 


0.01 


-0.06 




24 


3.27 


1.61 


-0.19 


0.13 




36 


8.20 


4:49 


-0.34 


0.45 




48 


21.39 


10.19 


-0.38 


1.74 


600 km N 


12 


1.94 


-0.40 


0.01 


-0.06 




24 


3.25 


1.61 


-0.19 


0.13 




36 


8.21 


4.51 


-0.34 


0.44 




48 


21.53 


10.28 


-0.40 


1.74 


400 km N 


12 


1.95 


-0.41 


0.02 


-0.06 




24 


3.21 


1.61 


. -0.18 


0.13 




36 


8.19 


4.51 


-0.34 


0.44 




48 


21.58 


10.33 


-0.41 


1.73 


200 km N 


12 


1.95 


-0.41 


0.02 


-0.06 




24 


3.17 


1.59 


-0.18 


0.13 




36 


8.13 


4.50 


-0.33 


0.44 




48 


21.53 


10.30 


-0.41 


1.72 


60 km S 


12 


1,94 


-0.42 


0.02 
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Figure 4.24 Lower perturbation amplitudes after I h and 4S h for the jet streak 
cross-track phase difference cases. 
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• Cross-track superposition is not extremely important in the jet streak case since 
the surface cyclogenesis occurs at nearly the same location relative to the jet 
streak regardless of the initial meridional location of the lower-level 
perturbation. The most favorable position for the initial low-level perturbation 
is at the location of the strongest low-level baroclinity. 

• Unlike previous experiments, the most favorable along-stream tilt for upper 
perturbation development is not the same as the most favorable tilt for the 
surface low. When a jet streak is present in the model, the upper perturbation 
grows fastest when the jet streak and the surface trough are nearly vertically 
stacked. 

E. SUMMARY 

A channel version of NORAPS that is frictionless-, adiabatic, quasi-linear and has 
a constant low-stability lapse rate (8 K km' 1 ) is used to investigate the superposition of 
perturbations centered at 400 mb and 800 mb. The simplified model is chosen so that 
different variables can be reliably interpreted. This approach has some limitations in 
direct application to real-world' cyclogenesis. Some additional variables that are not 
included in this research are discussed below. 

Static stability. The stability is set at a relatively low value to ensure that it is 
favorable for cvclogenesis. In actual cases, the low-level static stability may not be 
favorable, or the region in which it is favorable may be limited in horizontal or vertical 
extent. Additionally, the stability may become less favorable with time. A more stable 
atmosphere would reduce the growth rate for all cases and would result in a longer 
wavelength of maximum growth, but would not change the conclusions regarding 
superposition. The relative growth rates in the various cases would remain the same. 

Horizontal scale. A match in horizontal scales of the upper and lower-level 
features is assumed. Farrell (1985) used an initial value approach to show that jet-level 
and surface features that are matched in size exhibit significant growth when properly 
superposed. However, the horizontal scales of translating jet-level features and existing 
surface low pressure centers are often different. While disturbances with wavelengths 
as short as in these experiments are common in a low static stability maritime 
environment, a wide range of upper-level features is potentially available for coupling. 
If a good scale-match between the upper and lower perturbations is not achieved, 
growth will be reduced in all cases. Provided a reasonable match is made between the 
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scales at the two levels, the basic conclusions with respect to superposition are 
unchanged. 

Timing of initial superposition. These experiments involve the development of 
systems in which the upper and lower perturbations have been initially placed in 
certain superpositions. Many cases of maritime cyclogenesis involve the coupling of a 
translating upper feature with a quasi-stationary surface low that may be associated 
with a coastal front or a sea-surface temperature gradient region. This scenario would 
require a range of superpositions from unfavorable to favorable as the upper system 
moves toward the surface low. At the same time, both systems may be changing in 
size and intensity. Presumably, there is an optimum combination of timing and 
development that will minimize negative feedback during the unfavorable superposition 
period and maximize positive feedback when the systems are favorably coupled. 

Nonlinear diabatic processes. While the perturbations in these experiments are 

not infinitesimal, they are relatively small compared to the 'mean flow. Therefore, the 

growth rates and favorable superpositions indicated by these results are directly 

applicable only to. the early stages of cyclogenesis cases. As cyclones grow, they 
% * » 
interact with the mean flow, so the quasi-linear approach is not strictly valid for 

rapidly deepening cases. Also, diabatic processes are most important when 

perturbation amplitudes and associated vertical motions are. relatively large. Several 

experiments were repeated with physical parameterizations included in the model 

formulation. Although the growth rates of the perturbations were increased as 

physical processes became important late in the integrations, the relative sizes of the 

perturbations in the various cases remained the same. However, it is expected that the 

convergence; divergence pattern associated with the jet streak would be even more 

important when physical processes are included. In that case, the upper-level 

divergence in the right-entrance and left-exit quadrants would favor development below 

those areas. 

With these limitations in mind, these experiments suggest a number of 
conclusions that are applicable to actual cyclogenesis cases. 

• For a given mean flow, maximum growth rates are achieved when relatively 
large perturbations are present at upper and lower levels. When this is the case, 
the zonal superposition of the two perturbations is crucial. While favorable 
coupling results in the maximum growth rate possible, an unfavorable phase 
difference can cause decay of the features at both levels. 
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• The most favorable zonal tilt between the surface low and the upper-level wave 
is slightly less than a quarter wavelength. 

• The most favorable position of the surface low perpendicular to the jet stream 
depends on the characteristics of the mean flow. Several factors are important 
in the determination of the most favorable meridional tilt, including horizontal 
shear and low-level baroclinity. Increasing the horizontal shear moves the most 
favorable position northward while enhancing the surface temperature gradient 
moves it southward. The location of the most rapid growth early in a 
superposition episode is critically dependent on the location of the maximum 
surface temperature gradient. 
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V. SUPERPOSITION EFFECTS DURING GALE IOP-1 



In preceeding chapters, simplified models were used to investigate superposition 
between upper and lower-level disturbances in different mean flows. The conclusions 
from those experiments are consistent with general synoptic reasoning and can provide 
valuable guidance to operational forecasters and field experimenters. Cyclogenesis in 
nature can be thought of as a complex, nonlinear process involving interactions 
between time-dependent mean flows and finite-amplitude disturbances. Diabatic and 
adiabatic processes act on several temporal and spatial scales. Since the final 
development of the surface cyclone depends on the combined contributions of all these 
factors, it is difficult to isolate the effects of a single process. 

In this chapter, a quasi- Lagrangian diagnostic (QLD) approach is used to discern 
the various contributions in one case of cyclogenesis in the eastern United States. 
While a single case study is not sufficient to verify whether the superposition • 
conclusions can be directly transferred to operational use, a thorough diagnostic study 
of. one case can provide indications of the relative importance of superposition in 
cyclogenesis. The case chosen for this research occurred during the first intensive 
observational period (IOP) of the Genesis of Atlantic Lows Experiment (GALE) from 
19 to 21 January 1986. GALE was a large, cooperative experiment devoted to the 
study of the mesoscale and air-sea interaction processes involved in winter cvclogenesis 
over the east coast of the United States. The field portion of the experiment (15 
January to 15 March 1986) involved the collection of extensive observations on small 
time and space scales for an area centered on North Carolina. However, the data 
available for this research are only those observations in the normal synoptic network. 
Supplementary GALE observations are not included. 

IOP-1 was a case in which strong jet-level forcing existed, but the expected rapid 
surface cvclogenesis did not occur. Among the several possible reasons for the 
relatively weak surface cyclogenesis is the lack of favorable superposition between the 
upper-level disturbance and the incipient surface low(s). Relatively large static stability 
in the vicinity of the surface cvclogenesis may have also inhibited transfer of the 
available energy from the upper-level feature to the developing surface cyclone. The 
mass, vorticity, heat and moisture budgets for IOP-1 can be examined to describe the 
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reasons for the weak surface low development. It is desirable to also calculate the 
budgets based on model forecasts as well as those based on analyses because of the 
additional contributions to budget residuals that arise from inherent data uncertainties 
in the over-ocean observations and from the absence of direct observations of subgrid- 
scale processes. When forecasts are used, the budget residuals include only the effects 
of interpolation and truncation errors, in addition to any physical processes not 
included in the budget formulation. Additionally, forecast data are available as often 
as desired, rather than at just 12 h synoptic intervals. In this research, the forecast 
budgets are calculated every 6 h. 

If the budget based on the forecast is to be used with confidence, the model 
forecast must have accurately predicted the development of the system, and the 
common terms in the budget should be 'similar to those based on the analysis. In the 
following sections, the NORAPS operational forecast based on 1900 (for brevity all 
dates and times will be represented by four numbers - 1900 denotes 00 GMT 19 
January 1986) initial conditions will be verified. The forecast budgets will be compared 
to analysis budgets, and important features that provide clues to the dynamics of 
IOP-1 will be highlighted. Finally, the 1900 lower-tropospheric conditions will be 
altered to simulate conditions that would be more conducive to cyclogenesis based on 
the favorable superposition hypothesis. The NORAPS forecast from these simulated 
initial conditions has more rapid. lo\v-level cyclogenesis. Although this is only a single 
case, this simulation does suggest that absence of a favorable superposition could have 
contributed to the absence of a rapid cyclogenesis event in this case of strong upper- 
level forcing. 



A. NORAPS ANALYSES AND FORECASTS OF IOP-1 
1. Synoptic Discussion 

This section contains a brief description of the significant features of IOP-1. 
Pertle (1987) provides a more detailed discussion of the NORAPS analyses and 
forecasts. Although the data collection phase of IOP-1 began at 1812, the most 
significant development occurred during the 48 h period beginning at 1900, which is the 
interval studied in this chapter. 
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a. 00 GMT 19 January 1986 

The NORAPS forecasts used in this research are initialized from the 1900 
data base. At this time, a strong short wave trough is located over the Mississippi 
river valley from the Great Lakes to the Gulf of Mexico (Fig. 5.1a). A polar jet streak 
(60 m s* 1 maximum velocity) is near the upstream inflection point and a subtropical jet 
maximum (also 60 m s' 1 ) is east of SC. ' 

The short wave trough is clearly evident at 500 mb (Fig. 5.1b), where the 
cyclonic curvature of the wave combines with shear to create an elongated absolute 
vorticity maximum (24 x 10' 5 s' 1 ) over IA and MO. Positive vorticity advection (PVA) 
is occurring in and downstream of the short wave trough. In the southeastern U.S., 
southwesterly flow advects vorticity from a weak positive center over GA northward 
toward NC. 

The minimum sea-level pressure over the U.S. (1008 mb)- is over the Great 
Lakes, where a surface cyclone associated with the Midwest short wave is moving 
eastward (Fig. 5.1c). However, a frontal wave over KY (1009 mb) will soon become 
the dominant surface system, and is designated the "primary" low for this study. The 
lower trppospheric batoclinic region is in the cold air west and north of the frontal 
system. PVA aloft supports intensification of the primary wave. The ridge of high 
pressure over the Atlantic extends to Cape Hatteras, which results in southerly flow 
and warm advection .along the eastern seaboard, There is no indication at the surface 
of the 500 mb vorticity maximum over GA. 

A cross-section from northern WI across Cape Hatteras, NC to 32°N\ 
71°W (bold line in Fig. 5.1b) is plotted in Fig. 5.2. This cross-section is chosen to 
study conditions that affect the subsequent development of both the primary low and a 
secondary center that forms by 1912. The cross-section is approximately perpendicular 
to the mean tropospheric wind at 1900, so the perpendicular wind component is a good 
estimate of the actual wind. Two jet cores are present. The 40 m s' 1 wind maximum 
associated with the Midwest short wave on the polar front jet is at 325 mb over Lake 
Michigan. The other maximum (45 ms' 1 ) at 250 mb over the Gulf Stream is 
associated with the subtropical jet streak centered near 32°N\ 74°W. The winds 
seaward of Cape Hatteras between 200 and 300 mb do not appear to be balanced with 
the analyzed temperature field. Indeed, estimates of the maximum geostrophic wind 
speed based on the integrated thermal wind equation are only 30 to 35 m s' 1 . The 
apparent imbalance may be partly due to an ageostrophic component of the 
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Figure 5.1. NORAPS analysis at 00 GMT 19 January 19S6. (a) 250 mb heights (solid) 
in gpm and isotachs (dashed) in m s . (b) 500 mb heights (solid) in gpm and vorticitv 
(dashed) x 10' 5 s' 1 , (c) Sea-level pressure (solid) in mb. 1000-500 mb' thickness 
(dashed) in gpm and frontal positions from corresponding NMC analvsis. 
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Figure 5.2. Perpendicular wind component (solid) in m s" 1 and potential temperature 
(dashed) in degrees K in the 00 GMT 19 January 19S6 NO RAPS analysis along a 
vertical cross-section from WI to 3l°N, 72°W. 
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anticyclonically curved subtropical jet, and partly to analysis errors in the data-sparse 
region over the ocean. This analysis may be improved later when the dropsondes from 
the GALE data set are included. The low-level baroclinity, as represented by the 
potential temperature gradient and the vertical wind shear, is very w r eak everywhere 
east of the Appalachians. In particular, cold air damming and a low-level jet that were 
present in classic cyclogenesis cases such as the Presidents' Day cyclone are not present 
in this case. 

b. 12 GMT 19 January 1986 

After 12 h, the upper-level short wave has translated eastward and become 
sharper (Fig. 5.3a and b). The horizontal scale of the upper-level system (estimated 
from the distance between the trough and the downstream ridge) of approximately 
3000 km is similar to the systems in Chapters III and IV. Thus, the growth rate would 
be expected to -be large if the static stability is sufficiently low. The trough axis lies 
along 88° W with a 60 m s' 1 jet maximum upstream. The height gradient east of the 
trough is larger than it was at 1900. with 50 m s* 1 southwesterly winds over the 
Carolina coast. Another jet maximum of 60 m s’* is centered east- of the Delmarva 
peninsula. 

The 500 mb vorticity pattern reflects the sharper trough with an even- 
stronger positive vorticity maximum (28 x 10' 5 s’ 1 ) over IL (Fig. 5.3b). A weak 
vorticity center (14 x 10' 5 s' 1 ) persists just east of the Appalachians over the NC SC 
border, and weak PVA continues northeastward across VA. 

Both of these 500 mb vorticity maxima have associated surface features. 
The frontal wave has moved from KY to the OH/WV border and deepened to 1004 
mb. A second surface low pressure center that is associated with the coastal PVA is 
located in southern NC. Since this coastal system (called the "secondary" low in this 
study) is in a position that is theoretically ideal for a pre-existing surface cyclone, it 
should be the starting point for rapid intensification, provided other factors are 

favorable. However, this system forms too late (after the short wave has already 

induced the primary low), and is located too far to the southeast of the upper-level 

forcing for most rapid growth. Based on a 3000 km wavelength, the vertical tilt 

between the coastal low and the Midwest trough is approximately 105°. Furthermore, 
the secondary low is approximately 450 km equatorward of the axis of the strongest 
PVA associated with the Midwest trough. As a result, the sea-level pressure pattern is 
a broad, disorganized area of low pressure rather than a compact center. Southerly 
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Figure 5.3. As in Figure 5.1. except for 12 GMT 19 January 19S6. 
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winds ahead of the secondary low cause an onshore flow of warm, moist air from Cape 
Hatteras to Cape Cod, but the land surface is too warm to provide the contrast 
necessary for a significant coastal front, 'baroclinic zone. The strongest lower- 
troposphere baroclinic zones are west of the cold front associated with the primary low 
and much farther to the north over eastern Canada. 

The lower- tropospheric static stability is calculated for the layer between 
1000 and 500 mb using the formula from Holton (1979): 

<y = - (a/0) dQ/dp (5.1) 

The stability in the vicinity of the coastal low is relatively high (Fig. 5.4), since the 
1000-500 mb static stability for a standard atmosphere is approximately 10.6 x 10' 2 . 
(For comparison, the 1000-500 mb static stability of the low-stability atmosphere used 
in Chapters III and IV is 4.7 x 10‘ 2 .) The minimum stability in the NORAPS domain 
is less than 10 x 10‘ 2 cm 2 g’ 1 mb' 1 in the vicinity of the primary low’ cold front. The 
entire Atlantic coast region north of Charleston is more stable than the standard 
atmosphere. This stable environment would tend to reduce the feedback between the 
upper and lower features, and thus inhibit intensification; 
c. 00 GMT 20 January 1986 

The upper-level trough has continued to move eastward with little change 
in intensity (Fig. 5.5a). The 55 m s' 1 jet maximum in the downstream southerly flow. is 
centered over the XJ coast. At 500 mb. a closed low’ has developed over the southern 
Appalachians with a 28 x 10' 5 s' 1 vorticity maximum that is 150 km to the south of 
the height minimum. PVA continues from the GA coast to the Delaware Bay. The 
significant mid-troposphere cyclogenesis is evident in the GOES water vapor imagery 
(not shown) that has a symmetrical "doughnut" of relatively dry air around the 500 mb 
cut-off low. 

The strong cyclogenesis in the middle and upper troposphere is not evident 
at the surface (Fig. 5.5c). The surface centers are still separated by 450 km with the 
997 mb primary low in western PA and the 996 mb secondary low over the southern 
Chesapeake Bay. It appears that the diffuse lower tropospheric circulation has not 
been effectively organized into the thermal and vertical motion fields that can transfer 
significant energy into the low-level system. Additionally, little baroclinity is present in 
the vicinity of the surface lows. Strong onshore flow exists from XJ to Xova Scotia, 
but there is little warm advection because the low-level temperature gradient is weak. 
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Figure 5.4. Static stability in the 1000-500 mb layer for 12 GMT 19 January 1986. 
Contours are 2 x 10' 2 cm 3 g' 1 mb' 1 . 



d. 12 GMT 20 January 1986 

The upper cyclogenesis has continued (the 500 mb low has deepened 60 m 
in the previous 12 h) as it moved northeastward to the Chesapeake Bay (Fig. 5.6b). 
Because the 24 x 10' 5 s' 1 vorticitv center is nearly concentric with the height minimum, 
only small mid-tropospheric PVA remains. The primary’ surface low has been absorbed 
into the secondary system circulation, and the minimum sea-level pressure is 989 mb in 
northeastern PA. The secondary low is mature, and a third center is beginning to 
organize at the frontal triple point over RI. 

e. 00 GMT 21 January 1986 

By the end of IOP-1, the jet-level short wave, the 500 mb cut-off low and 
the surface cyclone have moved northeast to New England and have begun to fill as 
the system becomes vertically stacked from the surface to the tropopause (Fig. 5.7). 
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Figure 5.5. As in Figure 5.1, except for 00 GMT 20 January 19S6. 
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Figure 5.6. As in Figure 5.1. except for 12 GMT 20 January 19S6. 
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Figure 5.7 As in Fig. 5.1, except for 00 GMT 21 January 19S6. 
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2. NORAPS Forecast 

The NORAPS forecast from the 1900 analysis accurately predicted the 
movement, development and filling of the primary low during the first 24 h of the 
integration (Figs. 5.8 and 5.9a). It also predicted the development of the secondari- 
low when it was first identified as a low pressure center on NMC analysis (1918). The 
forecast of sea-level pressure for the secondary low was generally excellent (Fig. 5.9b). 
Deepening was slightly slow during the first 12 h, but the prediction of the minimum 
pressure and the trend of intensification followed by slow filling matches very well with 
the NORAPS and -NMC analyses. The forecast of the secondary low movement was 
reasonably accurate through the 36 h point in the integration (2012). During the final 
12 h of the prognosis, NORAPS forecast the surface low to become barotropically 
stacked beneath the 500 mb cutoff low and. decelerated the system at all levels. Thus, 
the position error is large at 2100. 

The 36 h forecast fields for 2012 are shown at three levels in Fig. 5.10a, b and 
c. The right panels (d, e and f) show the differences between the forecast fields and the 
NORAPS analyses valid at that time. The errors are relatively smdll and are due partly 
to errors in central height/ pressure and partly to errors in position. In most of the area 
of interest, the forecast cyclogenesis was too weak at 250 mb and 500 mb (Figs. 5.10d 
and e). The surface low is slightly weaker and slightly farther south than the analyzed 
center. For. the purposes of this study, NORAPS successfully forecast the development 
of a complex synoptic pattern with respect to timing, position and magnitude. If the 
common terms of the budgets of mass, vorticity, heat and moisture based on these 
forecasts also have a close resemblance to those in the budgets based on the analysis 
fields, then it can be assumed that the forecasts capture the basic physical processes 
that were important. In that case, some inferences can be made regarding the relative 
importance of other terms than can be only estimated from the forecast fields. 

B. QUASI-LAG RANGI AN DIAGNOSTIC BUDGETS 

The development of the primary and secondary IOP-1 lows is investigated using 
the technique of quasi- Lagrangian diagnostics (QLD), which was developed by 
Johnson and Downey (1975) and was subsequently used by Wash (1978), Calland 
(1983) and Rau (1986), among others. Budgets are calculated for a cylindrical volume 
that moves with the cyclone, which eliminates the effects due to translation. The 
remaining budget terms are related to the development of the system. 
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Figure 5.S. Analyzed (solid) and forecast (dotted) tracks of the two lows during lOP-1. 
Labels represent day and hour (i.e.. 1900 is 00 GMT 19 January 19S6). 
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Figure 5.9. Central pressure for the (a) primary and (b) secondary lows. NORAPS 
forecast (solid). NORAPS analysis (dashed) and N'MC analysis (dotted). X-axis labels 
represent day and hour (i.e., 1900 is 00 GMT 19 January 1986). 
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Figure 5.10. (a, b and c) As in Figure 5.1. except lor XORAPS 36 h forecast for 12 
GMT 20 January 1986. (d, e and 0 Height, pressure verification (forecast minus 
analysis) at (d) 250 mb, (e) 500 mb and (0 surface. Contours in (d) and (e) are 60 gpm 
and in (f) are 4 mb. 
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For this study, a budget radius of 4° lat. is chosen to capture the major part of 
the circulation involved in the primary and secondary lows and to limit the influence of 
one center on the other. At 2000, the centers are separated by approximately 4° lat. 
and the budget volumes have a significant overlap. Budgets of mass, vorticity and heat 
are calculated for both lows using NORAPS analyses and forecasts. Moisture budgets 
are presented for the forecast fields only as moisture is not analyzed in N'ORAPS. 
Although budgets are calculated separately for the primary low and the coastal low, 
the significant development in this case occurred in the secondary system. 
Additionally, the life-span of the primary low was only approximately 24 h. Therefore, 
the budgets for the coastal low are discussed in this section and the primary low 
budgets are presented in Appendix C. 

1. Mass Budget 

The mass budget is the simplest QLD budget because it includes no sources or 
sinks. In general, the mass budget for the coastal low has a secondary (vertical) 
circulation in which mass converges into the budget volume in the lower troposphere, 
moves upward in the volume, and then is exported by upper-level divergence. The 
circulation has a strong, but short-lived, maximum between 1921 and 2006 (Fig. 5.11a). 
Convergence is present in a very deep, layer (surface to 450 mb) and there is strong 
divergence at the jet level. The associated upward velocity maximum (not shown) is 
centered near the level of non-divergence (LND) at 2006 and occurs during the most 
rapid .deepening phase of IOP-1. However, this vertical circulation is not maintained, 
as the amplitude is diminished by 50 percent during the succeeding 12 h, and the LND 
moves toward lower levels. In the forecast, the convergence is somewhat stronger near 
the surface and the LND is not quite as high. However, the rapid increases and 
decreases in vertical circulation are similar to the analysis. 

The rapid decrease in the vertical circulation may be due to the diffuse nature 
of the surface pressure. The two lows cause low-level convergence over a wide area, 
but do not organize it into a well-defined column of upward motion. The jet-level 
divergence maximum at 2000 to 2006 decreases rapidly as the mass circulation 
diminishes. This decrease is associated with the movement of the off-shore jet streak 
away from the surface low by 2012 (Fig. 5.6a). If favorable superposition had occurred 
between the available upper-level forcing and a single surface system, the enhanced 
mass circulation probably would have been established early and maintained for a 
longer time. 
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Figure 5.1 1. Horizontal mass flux across the 4° lat. radius in the (a) analysis and (b) 
forecast of the secondary low. Solid (dashed) contours denote outflow from (inflow 
into) the 4° lat. radius budget volume. Contours are 5 x 10' 1 ' g s* 1 mb' 1 . 
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2. Vorticity Budget 

Development of the low can be described in terms of the change in absolute 
vorticity (£ a ) through the life-cycle (Fig. 5.12). The coastal low has relatively low 
values of at its formation, but it spins up rapidly by the end of the period. The 
most rapid increase occurs between 2000 and 2012 and is largest in the upper 
troposphere (around 350 mb). Although the vorticity is relatively constant in the 
vertical between the tropopause and 850 mb, the spin-up at the surface is only about 
half as strong as that aloft. If favorable superposition had been achieved between the 
PVA aloft and the surface center, the surface £ a might have increased nearly as fast as 
it did aloft. 

Examination of the terms in the QLD vorticity budget can yield clues to the 
processes that caused the spin-up in IOP-1. The vorticity equation can be written 

6$ a ;6t = - I/A j S a (V„- V„ ) dl - did p (a ^) _ _ (5.2) 

-? a (V • V) + k • (dvid p X 7 to) + Fr + R, 

where £ a is absolute vorticity, A is the area of a pressure surface within the budget 
volume, V n is the normal wind component, V Q is the cyclone's normal velocity 
component, to is vertical velocity and the overbar denotes an area-average. The term 
on the left side is the quasi-Lagrangian vorticity tendency. On the right side, the first- 
term is the horizontal transport (flux) of vorticity and the second term is the vertical 
flux convergence. • The vorticity divergence term (third term) is considered to be a 
source of positive vorticity. The fourth term is the tilting term and the fifth term is 
frictional dissipation. The last term on the right is the residual, which includes the 
sources and sinks not explicitly resolved by the other terms, plus the computational 
errors. If all the sources and sinks contributing to vorticity changes have been 
properly accounted for, and interpolation errors have not been introduced, the residual 
should be small compared to the leading terms. Interpolation and calculation errors 
are the major source of uncertainty in budget studies based on numerical model 
output. For the analyses, an uncertainty due to incomplete observations is also 
present. 

The leading budget terms (friction is small) are averaged through the lower 
troposphere (1000 to 500 mb) and combined in Fig. 5.13. The leading terms causing 
vorticity spin-up are the divergence term and the horizontal transport of C a . The 
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Figure 5.12. Area-averaged absolute vorticity over 4° lat. radius in the (a) analysis and 
(b) forecast of the secondary low. Contours are 2 x 10 * s . 
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Figure 5.13. Average over the lower troposphere (1000-500 mb) of the vorticity budget 
terms for the (a) analysis and (b) forecast of the secondary low. Time tendency (solid), 
lateral transport (dashed), vorticity divergence (dotted), vertical transport 
(dash-dotted>, tilting term (long-short dashed) and residual (dash-asterisk). 



118 



vertical transport and the residual term tend to decrease the lower tropospheric 
vorticity. 

The trends, magnitudes and relative importance of the various terms are quite 
similar between the analysis and forecast budgets. These vorticity budget calculations 
reinforce the earlier conclusion that NORAPS handled this case well. The residual in 
the analysis is approximately the size of the vertical flux convergence, but smaller than 
the horizontal transport and divergence terms. In the forecast budget, the residual is 
smaller than all three leading terms. Although smaller than the leading terms, the 
residuals are still relatively large, especially in the forecast budget. This error may be 
due to the parameterization of friction in the budget formulation. As this low occurred 
over land, it may be that friction is much larger than estimated in the budget. 
Consequently, a large negative residual is required to balance the vorticity budget. 

The divergence term is a source (sink) of positive vorticity below (above) the 
LND. For low'-level cyclogenesis, it is one of the leading terms in the spin-up in both 
the analysis and forecast (Fig. 5.14). Convergence in the lower levels concentrates the 
ambient C a over a deep column. A local maximum in the magnitude of the divergence 
term exists between 2003 and 2006. This maximum corresponds to the peak in the 
vertical circulation during that period. After 2009, the magnitude of the upper-level 
divergence decreases in both the analysis and forecast. The divergence term 
contribution to C a in the lower troposphere decreases in the analysis, but the spin-up in 
the lowest layers continues until 2015 in the forecast. If the superposition between the 
surface center and the upper-level forcing had been more favorable, it may be expected 
that the magnitude of the C a increases due to the divergence term would have been 
significantly greater. 

3. Heat Budget 

The heat budget during the life span of a cyclone provides direct indications of 
the physical properties that contribute to the development. The area-averaged 
advective form of the thermodynamic equation is 

ST/St = -1/A(\V-V 0 ) • V T - w^T/dp + <oa/c p + Q. (5.3) 

where T is temperature, a is specific volume, c p is the specific heat at constant pressure 
and an overbar denotes the area-average of a term. The term on the left is the quasi- 
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Figure 5.14. Vorticity divergence term within 4° lat. radius in the (a) analysis and (b) 
forecast of the secondary low. Solid (dashed) contours indicate vorticity increase 
(decrease). Contours are 5 x 10' 10 s'". 
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Lagrangian temperature tendency. The first two terms on the right side, are the 
horizontal and vertical temperature advections. The third term is the energy 
conversion from ZAPE to EAPE (CE in the energy budgets in Chapter IV). Diabatic 
heating (Q) cannot be directly calculated from the analysis fields, so it must be 
estimated from the budget residual (last term), which also contains the computational 
errors in the budget calculation. 

The energy conversion term and the vertical advection term have large 
magnitudes and opposite signs. For example, the energy conversion term reaches 
nearly 40° C d' 1 in the forecast. These two terms are often combined to form an 
adiabatic cooling term, <0(1^ - T), where Tj is the dry adiabatic lapse rate and T is 
the environmental lapse rate or vertical temperature gradient. - T) is constrained 
to be always positive in the model atmosphere, due to the dry convective adjustment 
scheme. As the mean vertical motion is upward in the region of the cyclone, i.e., 0) is 
less than zero, this term is negative and tends to cool the column. 

The adiabatic cooling term and the remaining three terms in (5.3) are 
vertically integrated from 1000 to 500 mb for each time period (Fig. 5.15). The time 
tendency of temperature based on the analysis is negative throughout the- period. That 
is, the cooling effects of the adiabatic term override the weak diabatic heating. The 
horizontal advection of temperature is small throughout the life-cycle, of the secondary 
low, which is partly due to the lack of warm advection associated with the onshore 
flow ahead of the system. If the cyclone had formed in a strongly baroclinic lower 
troposphere (hypothesized to be an important factor in rapid cyclogenesis), the 
horizontal thermal advection would probably have contributed to warming the 1000 - 
500 mb column. 

In the analysis, the dramatic increase in the vertical circulation that occurred 
between 2003 and 2009 is reflected in an increase in the magnitude of the adiabatic 
cooling, but no significant change in the diabatic heating. In the forecast, the 
magnitudes of both terms double between 1921 and 2003, then remain relatively large. 
The vertical velocity maximum in the forecast is slightly greater than in the analysis, 
which results in a 4°C d* 1 increase in the magnitude of the energy conversion term (not 
shown). Comparisons of the relative magnitudes of the heat budget terms between the 
analysis and forecast indicate that the diabatic heating and energy conversion terms are 
slightly greater in the forecast. However, the changes compensate each other and the 
time tendency of temperature is nearly the same in the forecast as in the analysis. The 
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Figure 5.15 Average over the lower troposphere (1000-500 mb) of the heat budget 
terms for the (a) analysis and (b) forecast of the secondary low. Time tendency (solid) 
horizontal advection (dashed), adiabatic cooling (dotted) and diabatic heatins ' 
(dot-dashed). 
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diabatic heating term, which is the primary source of heat for both lows, is investigated 
in more detail below. 

Since diabatic heating is determined as a budget residual, all calculation and 
interpolation errors are also included in this term. Based on XORAPS analyses, the 
maximum diabatic heating of the column containing the coastal low was centered near 
500 mb (Fig. 5.16a). A maximum in heating at this level does coincide with the 
maximum in vertical circulation peak value at 2006. The contribution to the lower- 
troposphere, column-average temperature tendency is partially countered by a 
maximum in diabatic cooling near the surface at the same time. The stratospheric 
diabatic heating maximum is attributed to computational errors. The budget based on 
the forecast of the secondary low is quite similar to the analysis (Fig. 5.16b). Diabatic 
heating is stronger (5°C d' 1 ) throughout the troposphere after 2003, but the general 
shape of the field is the same. The greater diabatic heating during the last half of the 
forecast is consistent with the differences found in previous budgets. Greater diabatic 
heating implies a larger vertical circulation, greater low-level convergence and a larger 
spin-up of in the lower levels in the forecast. The height of the diabatic heating 
maximum is lower (600 mb) in the forecast and the layer of diabatic cooling near the 
surface is absent in the model. 

The diabatic heating rate estimated from the heat budget residual for the 
forecast can be compared to the archived heating rates in the model output (Fig. 
5.17a). The model-predicted heating maximum. is slightly higher (500 mb versus 600 
mb) than the center diagnosed as a residual in the budget. The model predicts diabatic 
cooling below 900 mb, similar to the analysis budget. The differences between the 
budget-diagnosed diabatic heating and the model-predicted heating are generally small 
(Fig. 5.17b). Primary' areas of disagreement are near the surface, which arises because 
sensible heating is not included in the archived heating, and in the lower stratosphere, 
which is probably due to computational errors in the budget. 

In summary, the heat budget calculated from the NO RAPS forecast fields 
provides a fair approximation of the analysis budget. The relative contributions of the 
diabatic and adiabatic processes are similar in both budgets. Based on the verification 
of basic height fields and the comparison of mass, vorticity and heat budgets, it is 
concluded that the NORAPS forecast from 1900 provided a good description of the 
subsequent development of the important features in IOP-1. The moisture budget 
(following section) can only be made from the forecast fields, but it can be assumed 
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Figure 5.16. Diagnosed diabatic heating within 4° lat. radius in the (a) analysis and (b) 
forecast of the secondary low. Solid (dashed) contours indicate heating (cooling). 
Contours are 5°C d' 1 . 
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Figure 5.17. (a) Diabatic heating from archives of NORAPS during the secondary low 
forecast. Solid (dashed) contours indicate heating (cooling), (b) Difference between 
diagnosed and model-predicted diabatic heating. Solid (dashed) contours indicate 
diagnosed heating is larger (smaller). Contours are 5°C d' . 
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with reasonable confidence that it is an accurate representation of the moisture 
processes involved in the development of both lows. 

4. Moisture Budget 

The area-averaged moisture budget equation is similar to the heat budget, 

Sq/St = -1/A § q(V V o) dl * d ' d + E-P, (5.4) 

where q is the specific humidity, E is evaporation, P is precipitation and an overbar 
denotes the area-average of a term. The term on the left is the quasi-Lagrangian 
moisture tendency. On the right side, the first two terms are the horizontal and vertical 
moisture fluxes. Evaporation and precipitation are combined in the residual, or 
moisture source/sink term. The sign convention for the residual is negative (positive) if 
precipitation (evaporation) is occurring. This budget residual can be compared to the 
predicted moistening rate from the archived model output. 

As in the vorticity and heat budgets, the terms are integrated from 1000 to 
500 mb (Fig. 5..18a). The horizontal flux moistens the lower troposphere and the 
condensation/evaporation term acts as a sink (condensation dominates). The vertical 
flux of moisture also dries the lower troposphere, but the magnitude of this term is 
relatively small. The lower troposphere is becoming drier throughout the life-cycle of 
the coastal system because the combination of vertical flux and condensation exceeds 
the horizontal influx of moisture. The horizontal flux reaches a maximum at 2003, but 
quickly decreases to relatively small values during the next period. This indicates that 
the secondary low did not continue to tap into the moisture that presumably was 
present in the onshore flow northeast of the center. Although this may be due to the 
diffuse nature of the surface pressure distribution, the moisture from the Atlantic may 
have been too far northeast of the surface center, and did not provide an adequate 
source of latent heat to maintain a strong vertical circulation. 

The distribution of the moisture sink in the vertical (Fig. 5.18b) coincides with 
the location of the maximum heating in the forecast as diagnosed in the heat budget. 
The condensation maximum is centered near 600 mb at 2006. The model-predicted 
moistening rate (not shown) is quite similar to the moisture sink diagnosed from the 
moisture budget. There is no evidence of a stratospheric maximum as there was in the 
diabatic heating, which indicates that the "heating" at this level was probably- a result 
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Figure 5.18. (a) Average over the lower troposphere (1000-500 mb) of the moisture 
budget terms for the forecast of the secondary low. Time tendency (solid), moisture 
sink (dashed), horizontal flux (dotted) and vertical (lux (dot-dashed), (b) Moisture sink 
within 4° lat. radius for the secondary low forecast. Dashed contours indicate 
condensation in g kg' 1 d' 1 . 
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of errors in the heat budget, as previously stated. The net moisture sink throughout 
the lowest levels does not support a net cooling due to the evaporation of rain below 
850 mb that was indicated by the diabatic heating term in the analysis heat budget. 
However, the moisture sink distribution is consistent with the presence of diabatic 
heating only in the NORAPS forecast. In summary, the NORAPS forecast has more 
heating and moistening of the coastal low than is inferred from the analyses. It can 
not be determined whether this more vigorous circulation in the forecast is real or 
whether the weaker circulation inferred from the analysis is simply due to incomplete 
observations. 

C. SUMMARY AND HYPOTHESIS 

IOP-1 was a case in which strong upper-level forcing failed to cause significant 
low-level cyclogenesis. The short wave trough intensified rapidly at upper levels, which 
resulted in a very strong vorticity maximum (28 x 10* 5 s* 1 at 500 mb). However, the 
associated surface development was poorly organized. A low west of the Appalachians 
developed and filled within an 24 -h period as the center of action shifted southeast to a 
system that moved along the Atlantic seaboard. It appears that this coastal system 
was close .to the correct horizontal scale to provide an excellent initial circulation for 
rapid cyclogenesis, as in the classic Presidents' Day cyclone. However, the timing and 
superposition with the upper-level forcing did not appear to have been favorable, and 
the low-level environment was not conducive for intensification. Baroclinity along the 
coast was weak and the static stability was relatively high. 

There are many factors in this case that appear to have not been favorable for 
rapid cyclogenesis. In addition to improper superposition and a relatively stable low- 
level atmosphere, the location of the short wave over land may have detracted from the 
development of the surface cyclone. If the surface low had been located over the 
ocean, it would have been subjected to smaller surface friction and had a direct access 
to the energy available from the relatively warm water via fluxes of sensible and latent 
heat. From the point-of-view of the preceding chapters, the coastal low developed too 
late and too far southeast of the upper-level vorticity source. Therefore, the upper- 
level forcing induced development of a low to the northwest of the coastal low. When 
the coastal low and the upper feature did move to a position where mutual 
reinforcement was possible, the coastal system was too weak to achieve the maximum 
deepening that would have been possible if the lower-level system would have already 
had significant vorticity. 
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A more favorable superposition may have been achieved if the low-level 
baroclinity along the Atlantic coast had been larger on 19 January (i.e., a coastal front 
or cold-air damming against the eastern slope of the Appalachians). The presence of a 
strong thermal gradient would have encouraged more rapid development of the coastal 
low, which might have made it strong enough to better link with the shortwave aloft 
when it moved into a more favorable superposition around 00 GMT 20 January. 
Based on the conclusions in the "baroclinic" case in Chapter IV, a stronger low-level 
thermal gradient should shift the position of the most favorable surface development of 
the primary low equatorward of the jet core. The combination of a stronger coastal 
low and an equatorward shift in the most favorable surface position may result in a 
superposition that increases the surface development. Specifically, the energy 
conversion and diabatic heating terms in the heat budget should increase, along tvith 
the mass circulation. This hypothesis will be tested in the following section by 
artificially increasing the low-level baroclinity along the mid-Atlantic coast to provide 
modified initial conditions for NORAPS. If the superposition conclusions are correct, 
the primary low should shift to the southeast, the secondary low should dominate more 
quickly, and the deepening rate of the coastal cyclone should increase. 

D. MODIFIED FORECAST 

To test the above hypothesis, the initial conditions at 1900 are adjusted to make 
the lower troposphere less stable along the Atlantic seaboard. This is accomplished by 
inserting a low-level jet (LLJ) parallel to the coast form Charleston, SC to New York. 
To balance the changed winds, the mass field is altered in the lowest levels. While the 
forecast from the new initial conditions is quite similar to the operational NORAPS 
forecast, several significant changes support the hypothesis in the previous section. 

1. Numerical Procedure 

Analyzed winds between 700 and 1000 mb from Charleston to New York are 
replaced by a LLJ from the southwest. The jet width extends from the Appalachians 
to approximately 200 km seaward of the continent. A 300 km-wide "blend zone" on 
either side of the jet allows a smooth transition to the analyzed fields. Divergence is 
removed from the new wind pattern and the mass field is adjusted on pressure levels to 
balance the non-divergent winds. 

The operational NORAPS initialization scheme removes divergence, but does 
not balance the mass field on pressure surfaces as was done for the modified forecast. 
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To insure that this step does not introduce significant changes in the forecast, the 
operational forecast is re-created with mass balancing included. Forecast fields at all 
levels are nearly indistinguishable from the forecasts discussed in section B. The 
forecast minimum pressure (at 2018) is 1 mb higher in the balanced forecast and the 
maximum deepening rate is slightly smaller. Thus, any changes induced by the 
initialization scheme that includes mass balancing would tend to reduce the predicted 
intensification rate. 

2. Modified Initial Conditions 

Insertion of a southwesterly LLJ along the east coast and removal of 
divergence results in small changes in the winds elsewhere in the troposphere (Fig. 
5.19). The polar front and subtropical jet maxima are in the same locations, but are 
somewhat (3 to 7 m s' 1 ) weaker. Additionally, a local minimum near 500 mb over 
Wheeling, WV is more pronounced. The LLJ maximum is 32 m s' 1 at 620 mb over 
Cape Hatteras. The strong vertical wind shear below the LLJ indicates a lower 
Richardson number (R-) between Richmond and Cape Hatteras. To balance the 
increased winds, the boundary 7 layer temperature (not shown) is approximately 4°C 

, a 

cooler between Richmond and Cape Hatteras than in the analysis (Fig. 5.2). This 
weak cold dome simulates conditions that may result from moderate cold-air damming 
against the eastern slopes of the Appalachians. These revised initial conditions should 
be more favorable for development of the coastal low. Also, *the results from the 
"baroclinic" case in Chapter IV indicate that increasing the low-level shear moves the 
most favorable position of the lower perturbation toward the equatorvvard side of the 
jet. Thus, these changes in the initial conditions should move the primary low. which 
is induced by the vigorous upper-level forcing, toward the anticyclonic side of the 
upper-level jet. 

The reduced stability is evident in the lower-tropo spheric static stability at 
1900 (Fig. 5.20a). At this time, the coastal low is forming near the NC/SC border 
(asterisk in Fig. 5.20b). The stability in the modified initial conditions is lower than in 
the operational analysis throughout the eastern U.S. (Fig. 5.4). 

These modifications are expected to induce changes in the forecast that will 
lead to a more favorable superposition. Results from previous chapters indicate that 
superposition is most important when strong features exist at both levels. Although a 
large low-level vorticity center is not added in this case, more vigorous development of 
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Figure 5.19. As in Fig. 5.2. except for the modified initial conditions at 00 GMT 19 
January 1986. 
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F igure 5.20. Modified forecast for 12 GMT 19 January 19S6. (a) 1000-500 mb static 
stability (contours are 2 x 10'- cm 3 g* 1 mb' 1 ), (b) Sea level pressure (solid) in mb and 
lOOOoOO mb thickness (dashed) in gpm. 
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the coastal low is expected to create a surface center that interacts more strongly with 
upper-level forcing. 

3. Changes in the Forecast 

The influence of the modified initial conditions is evident at all levels and at all 
times of the forecast. In the lower troposphere, the position and relative intensification 
rates of the primary and secondary lows in the modified forecast change as expected. 
The track of the primary low is displaced to the southeast compared to the position in 
the operational forecast. At 1912, the lowest surface pressure is 165 km to the east- 
southeast (Fig. 5.20b). The lower stability and strong low-level shear also contribute 
to more rapid development of the coastal low. Cyclogenesis occurs approximately 6 h 
earlier, and the secondary low becomes the dominant system at 2000 (compared to 
2006 in the operational forecast). 

The cyclogenesis in the middle and upper troposphere in the modified forecast 
is slightly stronger than in the operational forecast. The 36 h modified forecast (valid 
at 2012) is shown in Fig. 5.21 for comparison with the analysis (Fig. 5.6) and 
operational forecast (Fig. 5.10). The 250 mb trough is approximately 80 meters deeper 
than in the operational forecast (Fig. 5.21d), and is actually closer to the analyzed 
height pattern at 2012. The modified 500 mb low (80 m deeper than the operational 
forecast) is also quite similar to the 2012 analysis. Since the upper-level trough is 
deeper, the height gradient is larger, and the jet-level wind maxima are approximately 8 
to 10 m s" 1 stronger than in the operational forecast. 

The maximum 12 h deepening of the coastal low is 14 mb in the modified 
forecast, which is larger than the rate in the operational run ( 1 1 mb) or in the forecast 
with balancing included (9 mb). . However, the rapid deepening in the modified forecast 
does not continue, and the minimum sea-level pressure of 989 mb is only 1 mb deeper 
than in the operational forecast (Fig. 5.21c). Thus, the changes induced in the low- 
level circulation are in the desired direction, but are not sufficient in magnitude to 
achieve the most favorable superposition between the coastal low and the strong 
upper-level positive vorticity. Despite the southeastward shift in the position of the 
primary low', the surface centers are still separated by more than 500 km. Thus, even 
after modification of the initial conditions, the along-stream and cross-stream distance 
between the upper-level vorticity center and the coastal low is too great to achieve the 
most rapid deepening. For an optimum superposition to occur in IOP-1, the entire 
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Figure 5.21. (a, b and c) As in Fig. 5.1. except for the modified 36 h forecast at 12 
GMT 20 January 19S6. (d, e and f) Difference between modified and operational 
forecasts (operational minus modified) at (d) 250 mb. (e) 500 mb and (f) surface. 
Contours in (d) and (e) are 20 gpm and in (f) are 2 mb. 
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middle and upper troposphere would have to be shifted eastward several hundred km 
with respect to the surface low. 

E. BUDGETS FOR THE MODIFIED FORECAST 

Quasi- Lagrangian diagnostic budgets are performed on the forecast from the 
modified 1900 conditions for comparison with the budgets based on the operational 
forecast. Since the modified forecast is quite similar to the original forecast, it is 
expected that the modified budgets will also be very much like the forecast budgets of 
the previous sections. 

1. Mass and Vorticity Budgets 

The mass circulation in the modified forecast of the coastal low (Fig. 5.22) is 

very similar to the operational forecast (Fig. 5.11b). The magnitudes of the low-level 

convergence and jet-level divergence maxima are nearly unchanged. The major change 

in the mass circulation' is that the strongest circulation occurs approximately 6 h earlier 

in the modified forecast. This earlier development in the circulation is consistent with 

the earlier formation of the coastal low (see Fig. 5.11a), and will be evident in each 

budget. The mass circulation in the modified forecast also decreases more rapidly at 

♦ 

the end of the forecast period. By 2018, the magnitude of the mass outflow (Fig. 5.22) 
is approximately 50 percent of that in the operational forecast (Fig. 5.11b). The low- 
level mass convergence is also smaller than in the operational forecast. By the end of 
the modified forecast, the mass circulation is minimal. Weak convergence remains 
below 800 mb, but divergence aloft is nearly zero. This rapid weakening of the vertical 
circulation occurs as the upper-level low moves directly over the surface center, which 
reflects the lack of favorable superposition between the levels. 

The earlier formation is also evident in the time-height cross-section of q a 
(Fig. 5.23). The vorticity in the modified forecast not only peaks earlier, but it also 
begins to decay by the end of the period. Although the upper-level vorticity evolution 
in the modified forecast is closer to the analysis (Fig. 5.12a), the desired conditions for 
a more favorable superposition between the upper-level and surface were not achieved. 
The earlier spinup and subsequent decrease in the cyclogenesis can also be seen in the 
lower-troposphere average vorticity budget (Fig. 5.24). The terms have the same 
general shape and magnitude as they had in the operational forecast budget, but the 
maximum values occur approximately 6 h earlier (as in the analysis) and the 
magnitudes decrease rapidly after 2015. The vorticity divergence and lateral transport 
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Figure 5.22. Horizontal mass flux across the 4° lat. radius in the modified forecast of 
the secondary low. Solid (dashed) contours denote outflow from (inflow into) the 4° 
lat. radius budget volume. Contours are 5 x 10' 13 g s' 1 mb* 1 . 

terms decrease so much that the quasi- Lagrangian time tendency of C a actually 
becomes negative after 2015 as the lower-level system begins to spin down. The spin- 
down occurs as the system becomes vertically stacked from the surface to the 
tropopause at 2018 (not shown). The decrease in £ a generation by the divergence term 
is consistent with the weaker mass -circulation at the end of the forecast. 

2. Heat and Moisture Budgets 

Indications of an earlier intensification and earlier decay are also found in the 
column-averaged budgets of heat and moisture for the coastal low (Fig. 5.25). The 
relative sizes and maximum magnitudes of all the budget terms are nearly the same as 
in the operational forecast (Figs. 5.15b and 5.18a), except that the maxima are 
achieved approximately 6 h earlier and the magnitudes are smaller by the end of the 
period. Although the amplitudes are larger and are achieved earlier, the levels of 
maximum heating and drying (Fig. 5.26) are nearly identical to the operational forecast 
(Figs. 5.16b and 5.18b). 
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Figure 5.23. Area-averaged absolute vorticitv over 4° lat. radius in the modified 
forecast of the secondary low. Contours are 2 x 10'- s' 1 . 



137 



20 




Figure 5.24. Average over the lower troposphere (1000-500 mb) of the vorticity budget 
terms for the modified forecast of the secondary low. Time tendency (solid), lateral 
transport (dashed), vorticity divergence (dotted), vertical transport (dash-dotted), tilting 
term (long-short dashed) and residual (dash-asterisk). 
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Figure 5.26. (a) Diabatic heating (°C d" 1 ) and (b) moisture sink (g kg" 1 d" 1 ) within 4° 
lat. radius in the modified forecast of the secondary low. 
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F. SUMMARY OF IOP-1 CASE STUDY 

IOP-1 was a case in which a vigorous upper-level PVA source was not favorably 
superposed with respect to the surface lows to create the conditions for rapid 
cyclogenesis. Additionally, the low-level environment was not favorable for 
cyclogenesis, as the stability was too high and low-level baroclinity was very weak. As 
a result, intensification of the surface low was minimal. The sea-level pressure pattern 
was broad and diffuse, with centers co-existing west of the Appalachians and along the 
Atlantic coast for much of the period. Although the coastal low formed in a position 
near the NC coast that is often favorable for subsequent rapid deepening, it was well to 
the east of the most favorable surface position with respect to the upper-level forcing. 
As in the cases discussed in Chapters III and IV, deepening does occur with these 
larger vertical tilts, but it is not as rapid. 

The NORAPS forecast from 00 GMT 19 January is quite credible through 36 h. 
The track, intensification and rapid filling of the primary' low' center west of the 
Appalachians are accurately predicted, and the genesis and deepening of the coastal 
low' are well-represented. This provides justification for also calculating the budgets for 
the forecast low's to compare wdth the budgets based on the analyses. 

Quasi- Lagrangian diagnostic budgets from the NORAPS analyses and forecast 
illustrated the lack of favorable coupling between the upper and low'er levels. The low' 
west of the Appalachians was induced by the strong upper-level processes, but 
developed slowiy, and failed to establish a significant vertical circulation before the 
upper center overtook the surface circulation and moved to the east. The resulting 
easterly tilt with height (unfavorable for intensification) caused rapid weakening of the 
primary low after 2000. 

The coastal low' formed to the southeast of the primary center as the primary low 
reached maximum intensity. As a result, the low'-level convergence and positive 
vorticity were dispersed over a wide area, and the coastal cyclone was still relatively 
weak as the upper-level center moved into a favorable superposition around 2000. 
Consequently, the maxima in the mass circulation, vorticity generation, diabatic 
heating and precipitation, w’hich were observed between 2000 and 2006, were not 
maintained, and intensification of the system ceased. 

Modifications to the initial conditions for the NORAPS forecast were made to 
simulate earlier development of the coastal low, and to move the most favorable 
position with respect to the upper forcing toward the coast. The modified initial 
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conditions simulated the cold-air damming and low-level jet that evidently contributed 
to the classical cyclogenesis in the Presidents' Day case. Although both low centers 
reacted as expected to the modified conditions, and the intensification rate of the 
coastal system was measurably increased, the modifications were not sufficient to bring 
about a truly favorable superposition or to sustain rapid deepening for more than a 
short time. Since creation of more favorable low-level environmental conditions along 
the Atlantic coast was not sufficient, more drastic modifications to reduce the vertical 
tilt between the upper and lower circulations would be necessary to achieve optimum 
conditions. 
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VI. SUMMARY AND CONCLUSIONS 



Rapid intensification of mid-latitude cyclones requires a synergistic interaction of 
factors on several scales and at all levels of the troposphere. A favorable superposition 
of PVA aloft over a pre-existing surface low facilitates this interaction. Whereas 
cyclogenesis can occur when small perturbations are present in a baroclinicallv 
unstable mean flow, it is hypothesized that rapid cyclogenesis begins with the 
superposition of relatively large perturbations. Therefore, this study investigates this 
favorable superposition between finite-amplitude perturbations in the upper and lower 
troposphere. The horizontal scales of the perturbations are the same at both levels and 
are equal to the most unstable wavelength of the baroclinicallv unstable mean flow. 

Experiments with idealized models (a two-level analytical model and an adiabatic, 
frictionless channel version of NORAPS) demonstrate that the initial alongstream 
phase difference between developing features is crucial when both the upper and lower- 
level perturbations are relatively strong. When either perturbation is significantly 
stronger than the other, the phase difference has less effect on subsequent deepening. 
The dominant perturbation in that case simply forces growth at the other level in a 
position that is favorable for intensification, and the (small) pre-existing perturbation 
at that level is of little consequence. However, the most rapid growth for given mean 
flow and static stability conditions occurs when relatively strong perturbations at both 
levels are favorably superposed. The most favorable phase difference (a westward tilt 
with height of slightly less than a quarter-wavelength) is consistent with synoptic 
experience. Westward tilts other than the most favorable also result in growth of the 
perturbations, but the growth rate is smaller. By contrast, when the initial tilt is 
eastward with height, the perturbations negatively interact and the amplitudes decrease 
at both levels until the perturbations shift to a westward-tilt orientation. 

The analytical model can only address the along-stream tilt between the 
perturbations. Thus, the simplified version of NORAPS is used to test variations in 
the north-south position of the surface low. The cross-stream tilts only modify the 
optimum growth rates that is determined by the zonal superposition. If the zonal tilt is 
near the optimum westward tilt with height (approximately 70° between 400 mb and 
the surface), the perturbations grow rapidly regardless of their initial placement north 
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or south of mid-channel. Each case has a most favorable cross-track position, which is 
determined by the mean flow conditions. If the low-level environment is highly 
baroclinic (vertical wind shear is large in the lowest layers), the most favorable pre- 
existing disturbance position shifts toward the south. If the horizontal wind shear in 
the middle and upper troposphere is enhanced (by increasing the zonal mean velocity 
of the jet stream), the most favorable disturbance position shifts toward the north. 
Regardless of the initial disturbance position, the developing surface center shifts in 
time toward the most favorable cross-track position as it intensifies. 

When the meridional tilt is varied, the relative growth rates of the surface 
disturbances change with time. The perturbation that is nearest the most favorable 
position for a given mean flow has the largest baroclinic energy conversions after 24 to 
36 h, but the growth rate early in the integration depends directly on the position of 
the disturbance relative to the zone of low-level surface baroclinitv. The nearer the 
perturbation is positioned to the largest surface temperature gradient, the larger the 
growth early in the integration. 

When the upper-level wave perturbation is replaced as the source of PVA aloft 
by a jet-streak, the conclusions with respect to superposition are unchanged. The 
growth rates of the surface perturbations depend directly on the initial 2onal phase 
difference between the streak ‘and the surface trough. As in the experiments with an 
upper-level wave, the perturbation that has the same meridional position as the 
strongest low-level baroclinic zone experiences the most rapid growth at the outset of 
the integration. However, the cross-stream position is less important when a jet streak 
provides strong upper-level forcing, because the streak effectively shifts the surface 
perturbation to a favorable position (under the left-exit quadrant) early in the 
integration, regardless of where it is initially positioned. 

Since these experiments are conducted with simplified (linear, adiabatic, 
frictionless) models, caution must be used in applying the results to the rapid 
cyclogenesis hypothesis. However, the results appear to be consistent with synoptic 
experience. 

An actual case of east coast cyclogenesis that occurred during IOP-1 of GALE 
(19-21 January 1986) is investigated with the conclusions from the two-level and 
NORAPS channel model studies in mind. Rapid intensification did not occur at the 
surface in this case, even though a dramatic upper-level cyclogenesis did develop. Two 
surface lows did develop: one west of the Appalachians associated with a strong short 
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wave and PVA at upper levels, and the second along the Carolina coast. Although the 
coastal low subsequently dominated, it did not establish a well-developed vertical link 
with the upper-level forcing, and intensification was relatively w r eak and short-lived. 
Although this absence of a rapid surface cyclogenesis may have been due to a 
combination of factors, the hypothesis tested here is that it was due to a lack of 
favorable superposition. 

By modifying the initial conditions in the IOP-1 case, the low'-level environment 
in the region of the Atlantic coast w'as made more favorable for cyclogenesis. The 
forecast from the modified initial conditions included earlier development of the coastal 
low and a shift of the Appalachian low r toward the anticyclonic side of the upper-level 
jet. Although the intensification rate of the coastal low was increased compared to the 
operational forecast, a truly favorable superposition with the upper-level system was 
not achieved. Although the separation between the upper-level vorticity center and the 
surface low in this case is too large to completely simulate favorable superposition, the 
modified forecast supports the conclusions from the simplified model experiments. The 
southeastward shift in the track of the Appalachian low. indicates that the most 
favorable position for a surface disturbance moves equatorward in response to the 
increased vertical wind shear in the lowest levels. By forcing earlier grow r th of the 
coastal low, the lower feature is stronger in the modified forecast when the upper-level 
forcing moves eastward into a slightly more favorable superposition. Consequently, 
the intensification rate in the modified forecast is significantly larger. 

This study isolates the effects of variations in the superposition between existing 
upper-level and lower-level disturbances. It demonstrates that when vigorous upper- 
level forcing interacts with a surface cyclone of the same horizontal scale, the relative 
position between the tw r o is critical to subsequent development at both levels. 
Therefore, rapid intensification not only depends on the strength of upper-level and 
surface features, and the favorability of the low-level environment, but also on the 
proper superposition of these prerequisites. 
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APPENDIX A 

TWO-LEVEL MODEL STANDARD ATMOSPHERE CASES 



I. BAROCLINIC INSTABILITY IN A STANDARD ATMOSPHERE 

The two-level model is tested in a standard atmosphere to insure the validity of 
the method used in the low stability cases of Chapter III. For a static stability which 
corresponds to a lapse rate of 6.25 K km' 1 and a mean flow shear of 20 ms' 1 , the 
growth rate curve is given in Fig. 3.2. The maximum growth rate (2.38 day' 1 ) occurs 
for a wavelength of 4000 km, consistent with expectations from baroclinic instability 
theory. The shomvave cutoff for this static stability is at 2591 km. Perturbations 
shorter than this will be neutral rather than unstable. 

The first standard atmosphere experiment (Case IS) is made with a wavelength of 
4000 km, an initial streamfunction perturbation at level 1 of 10 nrs' 1 , and a level 3 
perturbation of 1 nrs' 1 . The second, Case 2S, is made with the same, perturbation 
amplitudes, but with a wavelength of 2600 km. The amplitudes and phase 
relationships of the growing upper and low r er-level waves are shown in Fig. A.l. 
Positive phase differences indicate a wave which tilts westward with height. 

The Case IS upper-level wave grows from 10 m 2 s -1 to more than 30 nrs* 1 in 48, 
hours. The corresponding low r er perturbation, wiiich is initially an order of magnitude 
smaller, growls very rapidly and becomes nearly as large as the upper after 48 hours. 
The growth of the vertically-averaged perturbation from 5.5 nrs* 1 to 29.4 nrs' 1 
coincides with the expected growth rate. The perturbations, which were initially in 
phase, shift so that a westward tilt w r ith height (65°) is established. 

Case 2S uses a perturbation wavelength just slightly longer than the shomvave 
cutoff. The reduced instability is reflected in the lower growth rate, but the 800 mb 
perturbation again grows to the same order of magnitude as that at 400 mb in 48 
hours. As explained by Hoskins et al. (1985), the westward tilt for this shorter wave is 
not as great as in Case IS. A 55° tilt is established after 8 hours, but it quickly 
decreases to less than 30°. 
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CASES IS and 2S 

Perturbation Amplitude 





Figure A.l. (a) Amplitudes of 400 mb and 800 mb perturbations, and (b) phase 
differences through 48 h for Cases IS (solid) and 2S^ (dashed). Initial streamfunction 
perturbations in (a) are 10 nr s' 1 at level 1 and 1 mV 1 at level 3. 
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2. THE EFFECT OF VARIATIONS IN THE MEAN FLOW SHEAR 

Cases 3S and 4S investigate the importance of the shear in the mean flow. In 
both cases, the amplitude and wavelength of the initial perturbations and the static 
stability are identical to Case IS. In Case 3S the shear between levels 1 and 3 is 
10 ms' 1 , or half that of Case IS. Comparing Fig. A.2a to Fig. A. la, it can be easily 
seen that the amplitudes of both the upper and lower-level growing waves are much 
smaller when the mean shear is reduced. Conversely, if the shear is doubled (Case 4S), 
the waves grow so fast that by 48 hours the magnitudes of the perturbations are nearly 
six times as large as those found in Case IS. Indeed, when the shear is 40 ms' 1 , it 
takes only 24 hours to grow a wave which is as large as that in Case IS. 

The basic patterns of the growing wave are similar, regardless of the shear in the 
mean flo.w. The wave tilts westward with height ("against" the shear) and the phase 
separation between the 400 mb and 800 mb features is approximately 65° (Fig. A. 2b). 
The magnitude of the shear only affects the growth rate, as seen in equation (3.14). 
The (constant) shear in the mean flow betw'een level 3 and. level 1 is directly related to 
the meridional temperature gradient, dl!dy, by the thermal wind equation. A larger 
shear indicates a larger meridional temperature gradient, and thus a greater source of 
APE. Thus, doubling the shear will result in a two-fold increase in the growth rate. 

3. SHORT NEUTRAL WAVES 

The development of waves which are too short to be baroclinically unstable is 
also investigated using the tw T o-level model. Initial conditions are the same as in Case 
IS (static stability corresponding to a lapse rate of 6.25 K km' 1 , Uj = 10 m s' 1 , 
vj/j = 10 nrs' 1 , y-j = 1 irTs' 1 ). Perturbation wavelengths shorter than the shortwave 
cutoff (2591 km) are tested and 96 hour level 3 growth and phase differences are shown 
in Fig. A. 3. 

The neutral waves do not maintain a favorable vertical connection between the 
perturbations. The phase continually changes (Fig. A. 3b), and the w r aves grow when 
the tilt is westward with height. As soon as the phase difference changes sign, the 
amplitude of the w r ave begins to decay. The closer the wavelength is to the critical 
value, the longer the initial growth period (Fig. A. 3b) and the larger the maximum 
amplitude (Fig. A. 3a). For a 2000 km w r ave, an oscillation between growth and decay 
is established. 
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For waves which are just shorter than the shortwave cutofT, the period of growth 
can be quite long. The 2500 km wave does not shift to a decaying phase until 92 
hours. Naturally, linear model assumptions break down long before this. Still, the 
two-level model indicates that waves just shorter than the shortwave cutoff can achieve 
significant early growth. 
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CASES 3S and 4S 

Perturbation Amplitude 





Figure A.2. As in Figure A. I, except for Cases 3S (solid) and 4S (dashed). 
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SHORT NEUTRAL CASES 

Lower Perturbation Amplitude 





Figure A. 3. (a) Amplitudes of S00 mb perturbations, and (b) phase differences through 
96~h for wavelengths (L) less than the shortwave cutoff. L = 2000 km (dotted), 2200 
km (dot-dashed), 2400 km (dashed) and 2590 km (solid). 
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APPENDIX B 

NAVY OPERATIONAL REGIONAL ATMOSPHERIC PREDICTION 

SYSTEM-NORAPS 

1. MODEL CHARACTERISTICS 

NORAPS, which includes a regional analysis and forecast model, produces high 
spatial resolution (typically 80 km in the western Atlantic version), short-term (36-48 h) 
numerical forecasts over a limited domain. The principle advantage of using a regional 
model as opposed to a global or hemispheric model is the small spatial scale features 
that are resolved. This resolution is very important for simulations of maritime 
cyclones, which characteristically are smaller than continental systems and depend on 
interactions between fields that vary dramatically over short (meso-a scale) distances. 
XORAPS has the additional asset of flexibility, as the grid is globally relocatable, and 
the user may specify the forecast area dimensions and horizontal/vertical resolution. 
Three different projections (Mercator, Lambert conformal, or polar stereographic) are 
available to minimize distortion in the tropics, mid-latitudes or polar regions. An 
additional feature'is the "terrain enveloping" concept, in which topography is calculated 
at a high horizontal resolution to incorporate the effects of the sub-grid scale features 
into the topographic field. A thorough discussion of NORAPS is provided by Hodur 
(1982 , 1984). 

NORAPS includes analysis, initialization, forecast and output components. The 
analysis component begins with acquisition of different types of data (radiosonde, 
satellite, land and ship reports) and quality control checks to determine data validity. 
The next step in the analysis is the interpolation of the observations to the model grid. 
The data fields for the model are the u and v wind components, temperature, 
geopotential, specific humidity, surface pressure, sea-surface temperature and terrain 
height. A regional update cycle is used in which the 12 hour NORAPS forecast serves 
as the first guess. A successive corrections technique is used for the NORAPS 
objective analysis to improve the first-guess fields of wind, temperature, and 
geopotential. An exponential weighting function takes into account the distance of the 
observation from the grid point and the number of observations surrounding each grid 
point. The analyses are performed at eleven standard levels (1000, 925, 850, 700, 500, 
400, 300, 250, 200, 150 and 100 mb). 
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The initialization phase approximately balances the mass and wind fields and 
suppresses the growth of large amplitude inertial- gravity waves that would contaminate 
the forecast fields. The static initialization procedure uses diagnostic constraints to 
relate the wind and mass fields. 

The forecast component is the heart of the NORAPS model and requires the 
majority of the computation time. The model uses the flux form of the primitive 
equations on a staggered grid scheme C (Arakawa, 1977). This grid scheme has 
excellent geostrophic adjustment properties and group velocity characteristics (Haltiner 
and Williams, 1980). The vertical coordinate is sigma, which orients all coordinate 
surfaces parallel to the surface terrain. Thus, the vertical velocity is identically zero at 
the lower boundary, even in the vicinity of mountainous terrain. The vertical structure 
of the atmosphere is normally represented in 12 discrete layers, although as many as 20 
layers may be specified. All variables (u, v, q and T) except vertical velocity are carried 
at the middle of each layer. A split-explicit time integration scheme is used to permit 
larger time steps for the slower meteorological modes, while still being able to predict 
all the gravity modes. The size of the time step is governed by the computational 

i 1 ' * 

stability criterion for the horizontal resolution selected. Fourth-order advection'is used 
for the prediction equation set to reduce errors in phase speed. 

One-way interactive boundary’ conditions are used to specify the time-dependent 
lateral boundary conditions on the XORAPS grid from the coarser mesh NOGAPS 
predictions. The NOGAPS solution forces the fme-mesh model, without the fine grid 
affecting the coarse grid solution. Operational time constraints require that these 
boundary conditions be derived from the most recent forecast, rather than the forecast 
from the same time. A method developed by Perkey and Kreitzberg (1976) is used to 
spatially interpolate the solutions near the boundary’ of the finer mesh. The NOGAPS 
time tendencies are blended with the NORAPS time tendencies over a distance of five 
grid points to dampen spurious reflections due to the change in grid spacing at the 
regional model boundary’. 

The output phase of NORAPS prepares the forecast data for interpolation to 
standard pressure levels. Output fields include winds, temperature, specific humidity, 
surface pressure, relative humidity, absolute vorticity, divergence, surface sensible and 
latent heat flux, and precipitation. 
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2. MODEL PHYSICS 



The model physics are crucial for the simulation of maritime cyclones. NORAPS 
includes representations of the following physical processes: (i) dry convective 

adjustment; (ii) surface friction; (iii) cumulus parameterization; (iv) large-scale 
precipitation; and (v) radiative transfer processes. 

a. Planetary Boundary Layer 

The planetary boundary layer (PBL), which is confined to the bottom three 
layers of the model atmosphere, is well mixed in temperature, momentum and 
moisture. Interactions between the lower boundary and overlying atmospheric layer 
provide sources and sinks for momentum, heat and moisture. The effects of the PBL 
should be included in any numerical model to physically simulate maritime cyclogenesis 
on the time scales of more than a few hours (Anthes et al., 1983). 

The N’ORAPS PBL parameterization follows Deardorff (1972). After the layer 
mean values of V, 0 and q are known, a bulk Richardson number (Ri b ) is computed to 
determine the stability of the PBL, 

I 3 ' 

Ri b = gh(0 v - e v )/U m 2 , (B.l) 

m s 

where 0 y is the virtual potential temperature, subscript s denotes surface values, and 
subscript m denotes mean PBL values. For unstable conditions (Ri b <0), i.e., strong 
winds, daytime heating over land or strong mixing, a- predictive equation for the PBL 
height (h) proposed by Stull (1976) is used. The rate of change of the PBL height is 
related to the surface sensible heat flux, mean PBL wind speed, the large-scale vertical 
motion and cloud-induced subsidence. For stable or neutral conditions (Ri b >0), a 
predictive equation for the PBL height after Nieuwstadt and Tennekes (1981) is used. 
Transfer (drag) coefficients for heat (Cq) and friction (C u ) are computed from 
empirical formulas that include the stability dependence via the Ri^. Surface fluxes of 
moisture, heat and momentum are computed using the bulk aerodynamic formulas, 
which assume the transfer coefficients are functions of Ri b and h. Seasonally 
dependent climatological values of albedo, sea ice, ground wetness and surface 
roughness are specified. A predictive equation for the ground temperature after 
Blackadar (1979) is used to model the lower boundary condition for the temperature 
over land. The sea-surface temperatures are assumed to be constant over the forecast 
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period, which is valid for short-range forecasts. The PBL is constrained to extend 
through at least the bottom layer of the model (approximately 40 mb) to avoid 
extrapolation problems in determining mean PBL quantities. 

b. Cumulus Parameterization 

The NORAPS model uses a modified version of the Kuo (1965) cumulus 
parameterization scheme. This version links the convection to the PBL by requiring 
moisture convergence in the PBL. By contrast, the original Kuo version required net 
moisture convergence in the entire column before convection was initiated. The 
moisture convergence is 



M t = 1/g V • (q m *V m Xl-<T pbl ) + P s (w'q') s , ’ (B.5) 

where the first term on the right side is the vertically-integrated moisture convergence 
and the second term is the surface moisture flux. Convection is assumed to occur 
when M t > 4.0 x 10*^ gm m'^ s'^ and the equivalent potential temperature decreases 
with height (conditionally unstable) from the PBL to the first model layer above the 
PBL. The final constraint is that deep convection cqnnot occur if the lifting 
condensation level (LCL)-is above the PBL. -This scheme partitions the moisture 

transport into two fractions:, the first (bM t ) serves to moisten the environment to 

saturation conditions through the cloud layer, and’ the second ((l-b)M t ) condenses and 
falls instantaneously as rain. 

In addition to large-scale advection, temperature and moisture changes at any 
level are caused by convective clouds. It is assumed that the temperature of the cloud 
is warmer than the environment. Cloud production, which is the ratio of the water 
vapor available to the water vapor needed to form the cloud, is computed for each 

gridpoint. The fractional cloud area is used to adjust the layer mean temperatures and 

moistures to account for the presence of clouds. 

c. Large-Scale Precipitation 

Large-scale (non-convective) precipitation can occur when supersaturation is 
achieved at any level. The excess moisture is allowed to fall into the next layer and 
increase the moisture content of that layer, or continue to fall if that layer is already 
supersaturated. Precipitation occurs only when the air is saturated from the cloud to 
the ground. Physical processes evolve more slowly and are not as severely restricted by 
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time step limitations as are the dynamics. The precipitation routines are only called 
every eight time steps (once every 30 minutes) for computational efficiency. The 
heating and moistening rates are distributed evenly over subsequent time steps until the 
next call to these routines. 

d. Radiation 

The incorporation of solar radiation into numerical models is essential for 
prediction of surface temperatures and the cooling rates at cloud tops that may deepen 
cloud layers. The radiation parameterization in NORAPS follows Katayama (1974) for 
short-wave radiation and Sasamori (1968) for long-wave radiation. 
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APPENDIX C 

IOP-1 PRIMARY LOW BUDGETS 



The low pressure center which formed in KY at 1900 (denotes 00 GMT 19 
January 1986) during IOP-1 and filled in western PA after 2000 showed only minor 
development compared to a second low that formed on the NC coast at 1912 (Chapter 
V). Selected figures from the quasi-Lagrangian diagnostic budgets of the first low are 
presented for completeness in this appendix. 

1. MASS BUDGET 

In the mass divergence pattern for the IOP-1 primary low (Fig. C.l), 
convergence is limited to the lowest layers. In the analysis budget, the level of non- 
divergence (LND), which is initially at 600 mb, lowers to about 650 mb. and then rises 
as the primary low merges with the secondary system at 2000. Divergence aloft is 
weaker than the low-level -convergence, but extends through a- greater depth. The 
average vertical motion (not shown) in the budget volume is upward at all times, with 
the peak velocity at the LND (near 6S0 mb). There are only small changes with time 
in the mass circulation after the initial increase in boundary layer convergence. 

As expected from the synoptic discussion in Chapter V, the mass circulation in 
the forecast (Fig. C.lb) is quite similar to that in the analysis. Surface convergence 
and upper-level divergence are slightly stronger in the forecast, but- the general pattern 
is much the same. The LND remains below 600 mb throughout the period, but an 
area of convergence near 400 mb is present at 2000, as the primary system begins to be 
influenced by the secondary low. As in the analysis, there is little change in the pattern 
or- strength of the circulation from the beginning to the end of the budget. 

2. VORTICITY BUDGET 

Time-height cross-sections of show that the spin-up at all levels in this system 
was small. The analysis (Fig. C.2a) has a slight increase near the tropopause and an 
even smaller increase at lower levels. The forecast (Fig. C.2b) has only a small increase 
below 700 mb. In both budgets, a minimum in near 450 mb persists throughout the 
life of the low. 

The column-averaged vorticity budget for the primary low (Fig. C.3) has a 
similar pattern as that for the secondary low (Chapter V), but the magnitudes of all the 
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Figure C.l. Horizontal mass flux across the 4° lat. radius in the (a) analysis and (b) 
forecast of the primary low. Solid (dashed) contours denote outflow from (inflow into) 
the 4° lat. radius budget volume. Contours are 5 x 10' 13 g s' 1 mb' 1 . 
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Figure C.2. Area-averaged absolute vorticity over 4° lat. radius in the (a) analysis and 
(b) forecast of the primary low. Contours are 2 x 10' 5 s' 1 . 
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Figure C.3. Average over the lower troposphere (1000-500 mb) of the vorticity budget 
terms for the (a) analysis and (b) forecast of the primary low. Time tendency (solid), 
lateral transport (dashed), vorticity divergence (dotted), vertical transport 
(dash-dotted), tilting term (long-short dashed) and residual (dash-asterisk). 
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terms are smaller. Leading spin-up terms are the horizontal transport of and the 
divergence term, while the vertical flux convergence and the residual term tend to 
decrease the lower tropospheric vorticity in both the analysis and forecast. 

The vorticity generation due to the divergence term for the analysis and forecast 
of the primary low is presented in Fig. C.4. With the exception of a slight increase in 
the lowest-level spin-up, the analysis has only minor changes with time. In the forecast 
(Fig. C.4b), the upper-level decrease becomes smaller as the boundary layer spin-up 
continues. 

3. HEAT BUDGET 

The relative importance of the different terms in the column-averaged heat 
budgets for the analysis and forecast (Fig. C.5) are similar to those for the secondary 
low. The temperature tendency is small during all time periods in both the analysis 
and forecast, which is consistent with the slowly-changing nature of the low. Adiabatic 
cooling of about 7° C d* 1 in the analysis is balanced by small amounts of warm 
advection and diabatic heating. In the forecast budget, the magnitudes of the adiabatic 
cooling and diabatic heating terms are increased to near 10° C d* 1 . This adiabatic 
cooling reflects the stronger vertical circulation in the forecast that was identified in the 
mass budget. 

The residual (inferred diabatic heating) in the primary low is nearly constant in 
time (Fig. C.6). The level of maximum heating is near 650 mb in both the analysis and 
forecast. The secondary - maximum in the stratosphere, which is evident at 1906 in the 
analysis and at 1915 in the forecast, is probably attributable to accumulated errors. 
The diabatic heating in the forecast (Fig. C.6b) has the same pattern as in the 
analysis, except the magnitudes are consistently -5 0 C d' 1 larger. 

4. MOISTURE BUDGET 

The last budget for the IOP-1' lows is the moisture budget, which can only be 
calculated for the forecast fields. The column-averaged moisture budget has a pattern 
similar to that in the secondary - low (Chapter V), with the horizontal flux providing 
moisture to the lower troposphere, and the moisture sink (precipitation) term providing 
the greatest amount of drying (Fig. C.7a). The moisture tendency is positive in the 
early period, as the relatively large horizontal flux provides more moisture than is 
removed by condensation and the vertical flux. As the horizontal flux decreases 
toward the end of the budget calculation, the tendency becomes negative and the 
budget volume begins to dry. 
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As in the secondary low, the distribution of condensation in the vertical 
calculated as a budget residual coincides with the maximum drying as diagnosed in the 
forecast heat budget. The largest magnitude in the moisture sink term is centered near 
650 mb (Fig. C.7b). 
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Figure C.4. Vorticity divergence term within 4° lat. radius in the (a) analysis and (b) 
forecast of the primary low. Solid (dashed) contours indicate vorticitv increase 
(decrease). Contours are 5 x 10' 10 s' 2 . 
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1903 1909 1915 1921 



Figure C.5. Average over the lower troposphere (1000-500 mb) of the heat budget 
terms for the (a) analysis and (b) forecast of the primary low. Time tendency (solid), 
horizontal advection (dashed), adiabatic cooling (dotted) and diabatic heating 
(dot-dashed). 
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Figure C.6. Diagnosed diabatic heating within 4° lat. radius in the (a) analysis and (b) 
forecast of the primary low. Solid (dashed) contours indicate heating (cooling). 
Contours are 5°C d' 1 . 
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Figure C.7. (a) Average over the lower troposphere (1000-500 mb) of the moisture 
budget terms for the forecast of the primary low. Time tendency (solid), moisture sink 
(dashed), horizontal flux (dotted) and vertical flux (dot-dashed), (b) Moisture sink 
within 4° lat. radius for the primary low forecast. Dashed contours indicate 
condensation in g kg* 1 d* 1 . 
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